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ABSTRACT 


Whistlers  and  related  very  low  frequency  radio  signals  are  guided 
in  ducts  of  enhanced  or  reduced  ionization  along  the  geomagnetic  lines 
of  force  of  the  earth's  magnetosphere.  The  signals  convey  information 
about  the  distribution  of  particles  in  the  plasma  through  which  they  have 
propagated  and  about  the  occurrence  of  wave-particle  interactions  in  the 
magnetosphere.  Direction-finding  on  such  signals  will  aid  in  locating 
the  ducts  and  measuring  their  temporal  drifts,  thus  making  an  iiT5)ortant 
contribution  to.  studies  of  magnetospheric  convection. 

The  signals,  although  narrowband  in  nature,  exhibit  wide  frequency 
excursions  in  the  1 to  10  kilohertz  range.  An  innovative  technique  is 
presented  for  tracking  these  frequency  excursions  in  real  time  to  produce 
a filtered,  quasimonochromatic  version  of  the  signal.  The  voltages  in- 
duced by  the  incident  wave  on  two  orthogonal  loop  antennas  and  a vertical 
monopole  antenna  are  processed  by  this  method.  The  filtered  signals  are 
then  cross-multiplied  in  a manner  analogous  to  a Poynting  vector  calcula- 
tion to  obtain  a continuous  indication  of  the  wave  direction  of  arrival. 

The  design,  construction,  and  laboratory  testing  of  a prototype 
instrument  are  described.  Field  testing  of  the  instrument  was  performed 
at  Stanford,  California,  and  Roberval,  Quebec,  Canada.  Operating  at 
fixed  frequencies,  the  direction-finder  produced  accurate  results  on  VLF 
transmissions  in  the  10  to  20  kilohertz  range  from  NAA,  NP6,  GBR,  and 
the  Omega  stations.  Well  defined  bearing  indications  were  also  obtained 
for  spherics  at  frequencies  in  the  1 to  10  kilohertz  region.  Successful 
frequency-tracking  and  di Mictioh-findi ng  were  demonstrated  for  a wide 
variety  of  whistler-mode  signals,  tncluding  whistlers,  cho  dis- 
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in  the  2 to  7 kilohertz  range  stimulated  by  the  VLF  transmitter  at 
Siple  Station,  Antarctica. 

Direction  of  arrival  indications  on  whistlers  were  consistent  with 
the  duct  positions  predicted  by  frequency  profile  analysis.  Evidence 
was  also  presented  of  temporal  changes  in  the  direction  of  arrival  of 
signals  which  may  be  an.  indication  of  duct  drifts. 

This  instrument  provides,  for  the  first  time,  a continuous, 
immediately  readable  record  of  the  direction  of  arrive ' of  whistler-mode 
signals.  It  could  also  be  applied  to  the  automatic  acquisition  of 
whistler-mode  signals  and  measurement  of  their  (implitude  and  frequency 
variations. 
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A.  PURPOSE 

The  purpose  of  the  research  described  herein  is  to  develop  apparatus 
and  techniques  for  making  direction  of  arrival  measurements  on  whistlers 
and  other  very  low  frequency  (VLF)  signals  in  a field  environment. 

B.  BACKGROUND 

1.  Signals  of  Interest 

Whistlers  are  bursts  of  VLF  electromagnetic  energy  which  originate  in 
a lightning  stroke  and  are  guided  along  geomagnetic  lines  of  force  to  the 
opposite  hemisphere.  Their  character  is  that  of  a sweeping  oscillation 
descending  from  the  upper  audio  frequencies  < the  lower  audio  range  in  the 
course  of  a second  or  two;  the  frequency  variation  with  time  is  a result  of 
the  dispersive  propagation  characteristics  of  the  magretoplasma  through 
which  the  whistlers  propagate . It  is  generally  agreed  [Helliwell , 1965] 
that  the  whistlers  are  guided  in  "ducts'*  of  enhanced  or  reduced  ionization 
within  the  magnetosphere.  The  frequency -time  profile  of  a whistler  conveys 
considerable  information  about  particle  densities  in  those  regions  of  the 
magnetosphere  through  which  the  signal  has  passed.  As  a result , intensive 
scientific  investigation  has  been  underway  since  1951,  primarily  consisting 
of  the  recording  of  VLF  activity , spectrum  analysis  of  the  recorded  tapes , and 
interpretation  of  the  profiles. 
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While  the  frequency  pi’ofUe  analysis  can  pred.c^  the  equatorial  radius 
of  the  duct , and  hence  the  approximate  magnetic  latitude  of  its  endpoint 
j^Helliweil , 1965^  , no  longitudinal  infomation  has  been  avail^’e  other  than 
the  crude  inference  that  the  duct  exit  point  is  within  some  500  to  1000  km 
radius  of  the  recording  station . Studies  of  the  important  problem  of  temporal 
drift  of  the  ducts  [e.g. , Axford,  1969;  Carpenter  and  Seely,  1975 J , and 
experiments  to  detect  particles  precipitated  down  into  the  atmosphere  by  the 
passage  oi  waves  through  the  magnetosphere  require  a more  accurate  means  of 
locating  the  duct  exit  point . A single  direction-finder  determining  the 
azimuthal  angle  of  arrival  of  the  signal  could , in  conjunction  with  the  mag- 
netic latitude  analysis,  locate  the  ducts;  alternatively,  direction-finding 
measurements  at  multiple , spaced  stations  could  locate  the  ducts  by  triangula- 
tion, providing  further  data  on  the  accuracy  of  the  profile  analysis  as  a 
predictor  of  the  magnetic  latitude  of  the  duct  exit  point . These  opportunities 
provide  the  motivation  for  the  present  research . 

In  addition  to  whistlers . the  VLF  spectrum  exhibits  other  unusal  signals 
which  are  less  well-understood  than  whistlers  and  which  are  generally 
classified  as  **VLF  emissions" . Chorus  and  triggered  emissions  are  prime 
examples  of  these  signals . Of  particular  current  interest  are  the  artificially 
stimulated  emissions  (ASE)  produced  by  the  100  kW  VLF  transmitter  opei'ated 
by  Stanford  University  researdiers  at  Siple  Station , Antarctica  [Helliwell 
end  Katsufrakls , 1974]  . These  signals  present  an  opportunity  for  investigation 
beyond  the  realm  of  natural  whistlers;  the  stimulus  energy  can  be  controlled 
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to  enhance  the  ei^ecto  under  study  end  reduce  other  phenomena  which  mask 
them . It  is  therefore  desirable  that  a whistler-mode  direction-finder  be 
capable  of  operating  on  these  signals  aiso. 

2,  Classification  of  Dire^ion  Finding  Techniques . 

All  directioii-finders  are  apeciaiiged  antenna  and  receiver  systems 
which  determine  the  direction  of  arrival  of  an  incident  electromagnetic  wave . 
ClassiJlcaLicas  caii  be  made  according  to  which  wave  properties  are  measured 
and  according  to  how  these  measurements  are  processed  to  determine  the 
direction  of  arrival . 

Consider  a general  plane  wave  propagating  in  free  space  toward  a 
direction-finder  site  as  illustrated  in  Fig,  1.1.  The  azimuthal  angle  of’ 
arrival  is  denoted  by  0 and  the  elevation  angle  by  p . The  wave  polarization 
is  described  by  an  ellipse  having  an  axis  ratio  R and  a tilt  angle  a of  the 
major  axis . Any  arbitrary  polarization  cem  be  described  by  these  parameters; 
note  that  | R j 11,  and  that  R>0  corresponds  to  right-hand  polarization , while 
R<0  denotes  left-hand  polarization,  (These  geometrical  conventions  are 
similar  to  those  employed  by  Bullough  and  Sagredo  [l973j  and  are  convenient 
both  for  visualization  and  for  comparison  of  the  new  system  with  existing 
techniques.)  The  electric  and  magnetic  fields  of  the  wave  are  expressed  as: 

^ ^ icot  **  ^ 

E^e  (cos  « 3||  -sin  a ) + jEgf’’  (-sin  a -cos  a Sj^)  1 

a.l) 

H = Re  (cos  ck  + sin  a )4-j  ^ E J (-sin  o:  a +cos  a )J . 


O 

w 


Notes: 

1.  , ajj  , a^  are  the  uni.  vectors  in  the  wave  coordinate  system; 
is  the  direction  of  propagation , In  lies  in  the  plane  of 

incideuce , and  is  mutually  perpendicular  to  and  . 

2 . and  are  scalar  qaantities  indicating  the  electric  field 

along  the  major  and  -^^nor  axes  of  the  polarization  ellipse, 
respectively . 

3.  B =E2,/E^. 

4.  H is  positive  for  right-hand  polarization. 

5 . With  the  x and  y axes  aligned  shown , the  +z  axis  must  point 
downward  to  maintain  a Hght-handed  coordinate  s/stem . 

FIGUnE  1.1  A GENERAL  FLANE  WAVE  PROPAGATING  TOWARD  A DIRECTION- 
FINDER SITE. 
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In  accord  with  the  usual  convention,  we  will  drop  the  and  Re[J symbols. 
Also  applying  the  relation  , the  fields  can  be  expressed  as 

E = [(COS  Q!  a j - sina  ) + jR  (-sin  a - cos  a )] 

(1-2) 

H = ~E  [(cos  a a,  +sino!  a ) + jR  (-sinaa  + cos  a a,,)]  . 

^ ^ If  X ((  ■ 


Measurement  of  the  complete  vector  electric  mid  magnetic  fields  at  the 
direction-finder  site  is  sufficient  to  determine  the  direction  of  power  fiow 
from  Poynting’s  theorem.  In  theory , antennas  could  be  constructed  to  in- 
dependently sense  E , E , E , H , H , and  H and  circuits  could  be  employed 
X y 2 X y 2 

to  execute  the  multiplications  and  additions  necessarj^  to  compute  P^,  , and 

P , completely  defining  the  propagation  direction . Since  the  signal  may 
have  propagated  via  several  reflections  between  ♦he  earth  and  ionosphere, 
the  elevation  angle  of  arrival  is  of  limited  usefulness  in  locating  the  source . 
Hence , ground  based  systems  are  often  intentionally  limited  to  the  determin- 
ation of  azimuthal  direction  and  only  two  or  three  field  components  are 
measured . Direction-finders  operating  on  the  power-flow  measurement 
principle  can  utilize  antenna  systems  which  are  small  relative  to  the  wave- 
lengths involved  and  hence  may  be  termed  ’’small- aperture”  direction-finders . 

In  contrast , ’’large-aperture"  direction  finders  utilize  antenna  arrays 
which  are  appreciable  in  size  relative  to  the  signal  wavelengths . One  popular 
type  is  the  "interferometer"  direction-finder,  in  which  multiple,  spatially 
separated  antennas  are  employed.  Consider  a system  in  which  antennas  are 
located  at  the  origin  of  the  coordinate  system  of  Fig.  1.1  and  also  at  some 


radial  distance  along  the  x , y , and  z axes . Each  of  the  outlying  antennas 
will  sense  a wave  phase » relative  to  the  wave  phase  at  thw  origin , that 
corresponds  to  the  projection  of  the  wavenormal  vector  along  that  axis , 
thereby  defining  the  direction  of  propagation.  Again , most  practical  ground 
based  systems  are  limited  to  the  measurement  of  x and  y components,  and 
hence  only  the  azimuthal  arrival  angle  is  determined. 

The  accuracy  of  the  large-^erture  system  is  proportional  to  the 
spacings  used;  the  performance  of  the  two  types  of  direction-finders  in  the 
presence  of  instrumental  imperfections  and  noise  is  shown  in  Appendix  A to 
be  comparable  when  the  large-aperture  system  is  constructed  with  antenna 
spacings  of  V2v.  This  principle  is  responsible  for  the  practical  preference 
for  large*<eperture  systems  at  the  high-frequency  band  and  above,  where 
spacings  larger  than  this  are  easily  achieved,  and  for  the  predominance  of 
the  small- aperture  approach  in  the  VLF  range , where  such  spacings  require 
large  amounts  of  real  estaie  and  the  long  distance  transmission  of  phase  data. 

As  stated  earlier,  direction-finder}  can  Uso  be  classified  according  to 
how  the  wave  measurements  are  processed  to  determine  the  direction  of 
arrival,  regardless  of  which  wave  measurements  are  made. 

Consider  the  general  system  configurations  in  Fig.  1.2.  In  the  first 
system,  (a) , parallel  channels  filter  the  antenna  signals  and  then  a calcula- 
tion is  made  in  a non-linear  processor  to  determine  the  direction  of  arrival . 
Contrast  with  this  the  second  system,  G>)  * in  which  the  antenna  signals  are 
first  combined  in  some  time-varying  manner,  passed  through  a single 
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FIGURE  1 . 2 GENERAL  CONFIGURATIONS  FOR  PROCESSING  WAVE 

MEASUREMENTS  TO  DETERMINE  DIRECTION  OF  ARRIVAI . 
(a)  Multichannel,  parallel  processing,  (b)  Single-channel,  sequential 
processing. 
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filtering  channel  ^ and  then  processed  to  cellulate  the  direction  of  arrival . 
When  the  filtering  operation  requires  hviman  intervention , the  broadband 
unfiltered  signals  must  be  recorded . and  in  this  case  the  single  channel 
system  displays  the  advantage  of  requiring  only  one  recording  channel  as 
compared  to  the  requirement  for  matched  multiple  channels  in  the  parallel 
system.  However,  signals  of  brief  duration  and  signals  exhibiting  fading  or 
modulation  at  frequencies  comparable  to  the  time-varying  combiner  operating 
fiequency  are  prone  to  errors  or  complete  loss  of  directional  data  in  the  single 
channel  approach.  The  great  majority  of  work  in  VLF  direction-finding  has 
been  based  on  examination  by  an  analyst  of  the  entire  spectrum  to  locate  the 
signals  of  interest . Hence , the  single-channel  approach , embodying  the 
rotating  goniometer  as  the  time-varying  combiner,  has  predominated. 

3.  Prior  Efforts  in  VLF  Direction-Finding. 

The  earliest  interest  in  direction-finding  at  VLF  [ockenden , 1954  ; 
Adcock  and  Clarke , 194?]  was  apparently  motivated  by  a desire  to  locate 
the  sources  of  atmospherics  (generated  by  lightning  strokes)  for  meteor- 
ological purposes . In  the  early  systems , the  voltages  Induced  on  a pair  of 
crossed  loop  antennas  were  applied  through  dual  amplifiers  to  the  deflection 
plates  of  a cathode  ray  tube.  Direction  of  arrival  indication  was  purely 
visual  and  practically  instantaneous . In  the  classification  framework 
presented  above , these  systems  would  be  described  as  small-aperture 
direction-finders  using  multi-channel  parallel  signal  processing.  The 
polarization  errors  which  can  occur  in  such  crossed  loop  small-aperture 


systems  when  the  incident  wave  is  not  purely  vertically  polarized  were  dis- 
cussed by  Homer  [l954] . Using  continuous  wave  signals  from  station  GBR, 
Homer  [l95?]  also  demonstrated  that  while  daytime  polarization  errors  were 
low  (less  than  10°) , the  nighttime  values  ranged  up  to  40° , being  most 
prominent  at  distances  of  about  300  km  and  less  prevalent  at  shorter  and 
longer  distances . 

Specific  attacks  on  the  whistler  direction-finding  problem  began  in 
1955  with  Crary  [l96lj  , and  Watts  [l959j  who  used  a rotating  goniometer  to 
combine  the  outputs  of  the  crossed  loops  into  a composite  channel . These 
systems  followed  the  '’single-channel”  philosophy , with  the  goniometer  serving 
as  the  time-varying  combiner,  because  it  was  necessary  to  record  the  entire 
VLF  spectrum.  The  broadband  tape  recordings  were  spectrum  analyzed  onto 
paper  or  film,  and  then  visually  examined  to  locate  the  whistlers , determine 
their  ”nuU  pattern” , and  compare  the  pattern  with  some  reference  to  measure 
the  direction  of  arrival . Some  whistlers  exhibited  nulls  while  others  did  not . 
As  a possible  explanation , Watts  pointed  out  that  downcoming  circularly 
polarized  signals  would  not  display  good  nulls  on  such  a system.  Another 
justification  for  the  absence  of  nuUs  in  ‘‘swlshy”  whistlers  was  that  they  would 
probably  be  composed  of  many  elementary  whistlers  arriving  from  different 
directions . Ellis  [i960 j found  incomplete  nulls  when  measuring  VLF  noise 
at  S kHz  and  advanced  a similar  explanation:  the  source  must  be  of  wide 
angular  extent . The  most  recent  efforts  with  the  goniometer  and  crossed 
loops  have  been  published  by  Bullough  and  Sagredo  [l97sj  . Their  data  show 
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wide  variability  in  the  measured  direction  even  on  a single  whistler , but 
when  averaged  over  several  whistlers  judged  from  frequency  profile 

0 

analysis  to  have  the  same  path , the  RMS  deviation  is  on  the  order  of  tlO  to 
±20° . This  wide  variability  in  the  apparent  direction  was  predicted  in 
Crary's  [lOOil^  piupagation  analyses , along  with  indications  that  averaging 
the  data  over  frequ'.  icy  could  achieve  considerable  reductions  in  the  error . 

Tracing  the  history  of  *Targe-aperture"  systems  at  VLF , it  appears 
that  the  application  of  this  technique  for  whistlers  was  first  suggested  by 
Croiribie  [l955j  who  recommended  spacings  of  0.15 X , which,  interestingly , 
compares  closely  with  the  spacing  (X/2ir)  at  which  the  large- aperture  and 
small-aperture  techniques  display  comparable  accuracy  in  the  presence  of 
noise , The  first  report  on  an  actual  application  of  the  technique  is  given  by 
Delloue  [lOOOJ  , who  used  spaced  loops  and  a radio  relay  link  to  transmit  the 
phase  information . He  reported  that  the  direction  of  arrival  tended  to  "circle 
about  tho  magnetic  field  line"  and  that  the  polarization  varied  rapidly . 

Another  large  - ^erture  system  has  been  built  by  tlie  National  Bureau  of 
Standards  at  Brighton,  Colorado,  as  described  by  Idefley  et  al.  |^196lj  , with 
claimed  accuracy  of  under  1°;  however,  it  is  adapted  only  for  used  with 
atmospherics . The  antennas  are  placed  at  the  vertices  of  an  eqtiilateral 
triangle  with  sides  4 miles  in  length  and  the  signals  are  conveyed  to  the 
central  location  by  transmission  lines . 

Of  particular  interest  in  this  research  are  the  efforts  in  small-aperture , 
multi-channel  systems . Cousins  [l972j  proijosed  a method  for  combining  ihe 
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signals  from  crossed  loops  and  a vertical  antenna  which  would  yield  no 
polarization  error,  provided  the  antenna  responses  were  perfectly  in  agree- 
ment with  theory  and  the  incoming  signals  was  not  purely  vertically 
polarized!  He  recorded  the  three  antenna  signals  on  a multi-channel  tape 
recorder,  then  digitized  the  spectra,  and  performed  the  filtering  and  pro- 
cessing on  a computer  in  an  Interactive  mode  with  a human  operator.  Tsuruda 
and  Hayashi  [l974]  have  described  and  built  a system  using  a similar 
principle , but  employing  analog  processing  instead.  These  techniques  are 
characterized  by  the  need  for  multichannel  phase  and  amplitude  matched 
recording  capability;  human  intervention  is  required  to  locate  the  frequency- 
time profile  of  the  signal  and  a^ust  the  filtering  to  extract  it  before  the 
direction-finding  processing  is  performed. 

C.  APPROACH  SELECTED 

With  the  goal  of  this  research  being  the  development  of  apparatus 
suitable  for  field  usage , considerations  of  practicality  entered  the  design 
an  early  stage.  At  very  low  frequencies , the  small- aperture  antenna  system 
is  much  more  attractive  than  the  large- aperture  array  with  its  attendant  real 
estate  requirements  and  the  difficulty  involved  in  conveying  the  antenna 
signals  long  distances  to  the  processing  site . In  view  of  this  consideration , 
and  the  predominance  of  crossed  loop  antenna  systems  at  existing  VLF  re- 
ceiving sites,  the  small-aperture  approac’'  was  chosen.  The  more  fruitful 
area  for  improvement  was  found  to  b.  in  overcoming  the  limitation  of  existing 
t 

For  a discussion  of  this  point  see  Section  II . C . 5 . 
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VLF  data  gathering  and  signal  processing  techniques . Current  methods 
customarily  employ  either  broadband  recording  of  the  entire  spectrum  with 
subsequent  examination  of  spectrum  records  by  an  analyst , or  tixed  frequency 
narrowband  reception.  The  latter  method  can  capture  only  a fraction  of  the 
typical  whistler  or  artificially  stimulated  emission  in  the  course  of  its  wide 
frequency  variation  with  time.  A visible  gap  between  these  techniques 
exists:  broadband  analysis  provides  only  rough  indications  of  intensity  as  a 
function  of  frequency;  narrowband  techniques  produce  detailed  data  on  the 
signal  characteristics , but  capture  only  a fraction  of  the  desired  signal . If 
an  instrument  could  be  developed  which  would  automatically  detect  and 
filter  a frequency-varying  VUF  signal  from  the  background  energy  ^ real 
time,  this  gap  would  be  bridged.  The  following  benefits  would  be  expected: 

a . the  multichannel  processing  technique  of  Fig . 1 . 2 (a) 
would  become  feasibi'' , producing  a continuoua  record 
of  direction  of  arrival  uncompromised  by  fading  or  signal 
modulation,  and  U8id}le  for  signals  of  any  duration; 

b . a greater  amount  of  data  about  temporal  and  frequency- 
dependent  variations  in  the  apparent  direction  would  be 
produced , and  improved  accuracy  due  to  averaging  out  of 
errors  over  the  entire  signal  profile  would  be  expected; 

c.  automatic  extraction  of  other  signal  parameters  such  as 
foequenoy  and  amplitude  would  be  made  practical;  and 
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d . the  burden  on  the  human  observer  of  searching  and  scaling 
spectral  records  would  be  reduced. 

With  this  motivation , a new  approach  to  VLF  signal  data  gathering  was 
conceived , following  the  block  diagram  of  Fig.  1.3. 

D.  CONTRIBUTIONS  OF  THE  PRESENT  WORK 

The  most  important  contributions  of  this  research  are  summarized  as 
follows: 

1 . A new  and  innovative  technique  for  filtering  and  extracting 
direction  of  arrival  information  from  whistlers  and  other  VLF 
signals  is  proposed;  the  method  is  analyzed  and  compared  to 
alternative  approaches  in  theory;  a prototype  instrument 
embodying  the  new  technique  is  desigried  and  constructed 
and  deployed  in  an  actual  field  environment . The  practicality , 
efficiency , and  utility  of  the  approach  under  these  conditions 
are  decwinstrated . 

2 . The  direction-finding  data  produced  are  of  hitter  accuracy  than 
that  fiom  existing  approaches . displaying  less  polarization  error 
and  no  mnbigulty;  the  temporal  resolution  is  continuous  as 
compared  to  10  or  25  sample-per-seoend  limitations  of  the 
conventional  goniometer.  These  advantages  permit  more  detailed 
direction  arrival  studies  of  the  substructure  within  whistler 
events , including  precursors  and  triggered  emissions , as 
demonstrated  in  the  observed  data. 
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FIGURE  1.3  NEW  APPROACH  TO  VLF  DIRECTION  FINDING  AND  DATA  GATHERING 


3 . Examples  o*  temporal  changes  in  the  direction  of  arrival 
of  whistler-mcde  signals  are  found  which  may  indicate 
drifting  of  msgnetospheric  ducts . 

4 . The  firs',  -iirection-finding  on  signals  received  via  magneto- 
spheric  amplification  from  Siple  Station,  Antarctica,  is 
demonstrated.  Duct  exit  points  indicated  by  the  direction-finder 
agfrea  with  those  of  whistlers  received  in  the  same  time  period. 
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II.  THEORY  AND  ANALYSIS 


A . FREQUENCY  -TRACKING  OF  VLF  SIGNALS 

Th'a  concept  of  adjusting  the  frequency  of  a receiver  or  filter  to  ficiiow 
the  frequency  of  a slgsxal  has  been  applied  in  other  contexts;  in  particular, 
the  technology  of  phase-locked  receivers  for  satellite  signsil  reception  is 
highly  advanced , as  described  by  Gardner  [l966 j . All  of  these  applications , 
however , involve  tracking  over  a narrow  frequency  range  relative  to  the 
signal  frequency , whereas  a whistler  at  VLF  may  sweep  over  a decade  range 
in  frequency . To  accommodate  such  an  excursion , new  concepts  in  the 
design  of  the  frequency-tracking  circuitry  must  b«'  introduced. 

1.  Tracking  Loop  Characteristics  and  Limitations . 

Two  conflgurations  for  achieving  automatic  tuning  of  a receiver  or 
filter,  open -loop  and  closed-ioop,  are  illustrated  m Fig.  2.1.  V/hen  the 
purpose  of  the  filtering  is  to  extract  the  signal  from  background  noise  and 
interference , the  open-loop  approach  is  clearly  unsuitable  as  the  frequency 
information  would  contain  these  undesirable  elements  and  tracking  would  be 
poor;  open-loop  operation  is  also  more  sensitive  to  small  drifts  in  circuit 
parameters . Thus , the  logical  choice  is  the  closed-loop  configuration , in 
which  the  frequency  information  is  extracted  from  the  filtered  signal  and  fed 
back  to  control  the  tuning  of  the  receiver.  Furthermore , the  most  practical 
way  to  achieve  electronic  tuning  of  multiple  filtering  channels  is  to  use  fixed 
filters  and  to  instead  translate  the  frequency  of  the  incoming  signals  by  an 
electronically  variable  amount , using  a voltage-controlled  local  oscillator 
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and  mixers.  Implementation  of  this  technique  is  illustrated  in  Fig.  2.2. 

A rigorous  analysis  of  the  operation  of  such  circuitry  is  beyond  the 
realm  of  linear  circuit  theory , because  of  the  presence  of  nonlinear  elements: 
mixers , voltage-controlled  oscillators , limiters  and  frequency  discriminators . 
However , by  proper  definition  of  the  system  variables  and  by  restriction  of 
the  range  of  operation,  a linear  model  can  be  produced.  The  voltage- 
controlled  oscillator  and  frequency  discriminator  can  be  modeled  as  linear, 
memoryless  voltage-to-frequency  and  frequency -to- volt  age  transducers . The 
mixer , which  generates  sum  and  difference  frequencies  from  its  input  signals , 
can  be  represented  as  a frequency  subtractor.  (The  sum  frequency  component 
is  not  passed  by  the  IF  filter.)  For  the  IF  filter,  the  problem  lies  in  deter- 
mining the  response  of  the  output  frequency  to  changes  in  the  input  frequency . 
Panter  [l965]  and  Baghdady  [l96l]  have  performed  analyses  demonstrating 
that  an  appropriate  linear  model  for  the  "modulation  transfer  function"  is  pro- 
vided by  mathematically  translating  the  IF  bandpass  filter  function  into  a 
low-pass  filter  function  centered  at  zero  frequency . With  the  benefit  of  these 
approximations , an  analog  to  the  system  can  be  constructed  as  in  Fig.  2.3. 

Changes  in  the  tuning  frequency  of  the  receiver  are  brought  about 
by  the  detected  frequency  error  Af  in  the  IF  passband.  The  error  which  must 
be  tolerated  to  produce  the  necessary  changes  depends  on  the  loop  feedback 
filter  used.  A pure,  memory  less  amplifier  for  feedback  requires  an  error 
proportional  to  the  tuning  offset  needed  (that  is , the  difference  between  the 
signal  frequency  and  the  receiver  quiescent  tuning  frequency) . With  a 
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FIGURE  2 . 3 LINEAR  MODEL  OF  FREQUENCY-TRACKING  lOOP . 


single  integrator  in  the  loop , the  error  is  zero  fer  any  constant  frequency 

sign^s  and  pn>portional  to  df/dt  for  changing;  frequencies.  With  double 

2 3 

integration , the  error  is  propo^ional  to  d f/dt  , and  so  on;  however,  the 
higher  the  order  pf  integration , the  slower  the  initial  response  of  the  loop 
to  the  transients  involved  in  acquiring  a new  signal . For  whistlers , a single 
integrator  cq>pears  to  be  the  optimum  choice . 

Let  us  now  address  the  question  of  loop  gain  and  stability . The  DC 
gain  constants  around  the  loop  may  be  lumped  into  a single  factor  K and 
incorporated  into  the  integrator  gain , giving 

Hj^(s)  = K/s  (2.1) 

as  the  LaPlace  transform  transfer  function  of  the  loop  feedback  filter.  The 
modulation  transfer  huictiou  for  the  IF  filter  depends  on  the  exact  shape  of  the 
filter  response;  however,  a good  approximation  for  most  practical  circuits  is 
given  by 


Hj(,)  = -2 


8-  * (BW^,^/2)8  . (BW^^^/2)- 


(2.2) 


Where  is  the  bandwidth  of  the  IF  bandpass  filter  in  radians  per  second. 

This  function  exhibits  a conjugate  pole  pair  at  the  natural  frequency  , 

with  damping  ratio  i . The  response  of  the  entire  loop  is  then  described  by 
the  roots  of  the  characteristic  polynomial, 
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s®  + (BW^^^/2)  8^  + (BW^^^/2)^.s  + ^ * (2.3) 

Applying  the  Routh-Hurwitz  stability  criterion  [fiSelsaand  Schultz,  1S69]  , the 
loop  response  will  be  stable  if 

K < • (2.4) 

Allowing  the  customary  stability  margin,  a convenient  value,  for  the  gain  is 

K * (2.5) 

(Hereafter,  BW.„  will  be  taken  to  represent  BW  .) 

Ir  HZ 

In  order  to  simplify  the  calculation  of  system  response  to  various  inputs , 
it  can  be  observed  that  the  loop  behavior  is  dominated  by  the  integrator  pole 
at  zero  frequency;  thus , although  the  complete  loop  transfer  function  must  be 
considered  in  determining  stability , a modified  loop  model  as  shown  in 
Fig.  2.4  is  adequate  for  determining  system  response.  Using  K=BWjp,  the 
differential  equation  fer  ihe  frequency  error  Af  is  given  by 

df/dt  - dAf/dt 

At  - gj;,-  (2.6) 

from  which  we  can  conclude  that: 

a.  for  a constant  frequency  input , the  steady  state  error  is  zero; 

b . for  a small  step  change  in  frequency , the  transient  response 
has  a time  constant  of  (l/BW^j,); 
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FIGURE  2 . 4 SIMFLIFEEU  SYSTEM  RESPONSE  MODEL . 
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c . for  a frequency  ramp  input , the  steady  state  error  is  eonsteoit 
and  equal  to  7™^.;  and 

d . the  maximum  rate  of  change  in  frequency,  or  fastest  ramp  v;hich 

can  be  tracked,  is  that  for  which  1 Af|=  (BW_/2) , in  order  that 

Ir 

the  signal  remain  within  the  IF  passband.  Therefore, 


Idf/dtl 


max 


BWjy/2. 


(^.7) 


2 . Tracking  in  the  Presence  of  Noise  and  Interference . 

Let  us  postulate  the  addition  of  white  noise  to  the  signal . The  noise 
power  spectral  density  shall  be  denoted  by  , and  the  signal  power  by  S . 
Then  the  signal-to-noise  ratio  in  the  IF  is , 


' -N  lw„  - ■ «.8) 

0 IF 

Gardner  |l966j  has  shown  that  passage  of  the  signal  and  noise  through  an 
amplitude  limiter  has  a moderate  enhancement  effect  on  signal-to-noise  ratio 
as  given  b y 


(S/N) 


LIM 


(S/N)jp 


1 + 2(S/N)jp 

■4/ir+  (S/N)jp 


Which,  is  well  approximated  by 


(2.9) 


(S/N)j^  = 2(S/N)jp 


(2.10) 
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for  all  (S/N)jp  greater  than  0 dB.  Applying  the  limited  signal  to  the  dis- 
criminator, Banter  [l96sj  has  shown  an  analysis  demonstrating  that  the  voltage 
output  of  the  discriminator  will  contain  a mean-square  fluctuating  component 


of 


e 


n 


1.64 


(S/N) 


Bw: 


LIM 


IF 


(2.11) 


which , combining  the  limiter  and  discriminator , gives 
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e 


n 


0.82 

(S/N)jj. 


(2.12) 


The  noise  voltage  may  be  equivalently  expressed  as  representing  a frequency 

error  due  to  noise , Af  . Since 

n 

the  RMS  value  for  Af^  is  given  by , 


(Af^) 

" RMS 


0.91 


BW 


IF 


(2.13) 


Thus , the  frequency  error  is  a aero  mean , fluctuating  component  with  an  RMS 
value  approximately  equal  to  the  bandwidth  divided  by  the  square  root  of  the 
signal-to-nolse  ri^io . 

At  moderate  to  high  signal-to-noise  ratios , the  presence  of  such  a 
noise  component  does  not  disrupt  the  basic  tracking  behavior  of  the  receiver. 
At  low  signal-to-noise  ratios . however . the  fluctuation  noise  exceeds  the 
range  of  the  linear  approidmationa  made  in  studying  the  loop  behavior.  A 
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rigorous  solution  to  the  behavior  of  the  loop  under  such  conditions  has  not 
been  demonstrated:  however,  a number  of  workers  have  applied  expvorimental, 
simulative . and  analytical  techniques  with  the  conclusion  that  a sharp 
threshold  exists  near  unity  signal-to-noise  ratio  below  which  tracking  is  not 
possible  [Panter,  1965]  . 

Let  us  now  study  instead  the  addition  of  a second . coherent  signal 
within  the  IF  passband.  This  interfering  signal  mig^t  be  due  to  another 
source  or  it  could  be  a component  fx>om  the  same  origin  which  has  propagated 
via  a different  path . The  frequencies  of  the  desired  and  undesired  signals  are 
denoted  by  f.  and  f„ , respectively , and  the  power  ratio  as  (S/I)._.  Both 
signals  are  passed  by  the  IF  filter  and  applied  to  the  limiter- discriminator. 

The  effect  of  the  limiter  in  this  case  has  been  shown  by  Baghdady  [^196l]  to 
be  a moderate  suppression  of  the  weaker  signal  in  the  presence  of  the 
stronger  one , approximated  by 


(S/I) 


LIM 


(S/I) 


IF 


1 

4(S/I)jp 


2 


for  (S/I),^j  >2,  which  in  turn  can  be  approximated  as 


(2.14) 


(2.15) 


Applying  this  factor,  the  ampUtude  of  the  interfering  signal  relative  to  a 
unit  amplitude  desired  signal  is 
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at  the  input  to  the  discriminator.  The  veotoral  combination  of  these  two 
signals  is  illustrated  in  Fig.  2.5.  The  phase  of  the  resultant  of  the  two 
signals  is  given  by 

'Aj  sln2ir  (fg -f^)  t ' 

.1  + Aj  cos  2ir  (fg  J • 


= 2ir  f,  t + tan 


-1 


The  interference  term , represented  by  the  second  part  of  the  expression  for 
, can  be  expressed  alternatively  as  a aeries , giving 

^ (-1)”'^^  (Aj)*' (l/n)  Sin  2im(f2-fj)t. 
tr=l 

and . taking  deriv^ves  to  convert  phase  into  instantaneous  frequency , we 
have 

^ (-1)°^^  (A  )“c08  2im(f  -f  )t  . (2.17) 

iv=l 

Thus , the  output  of  the  frequency  discriminator  will  contain  a steady 
component  correctly  representing  the  frequency  of  the  strong  signal  , plus 
a time-periodic  error  component  tif^  due  to  the  presence  of  the  weaker  inter- 
fering signal.  The  character  of  is  illuatrated  In  Fig.  2.6  [panter,  1965]  . 
From  the  form  of  the  expression  for  Af^  > a summation  of  harmonic  cosinusoids 
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FIGURE  2.5 
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7ECT0RAL  COMBINATION  OF  A DESIRED,  UNIT-AAi?UTUDE 
SIGNAL  AT  FREQUENCY  AND  AN  INTERFERING  SIGNAL  OF 

AT  FREQUENCY  • 
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of  fundamental  frequency  conclude  that  the  time-averaged 

frequency  error , Af. , is  zero.  Furthermore , for  less  than  i , the  RMS 
frequency  error, is  well  approximated  by  . 

For  values  of  less  than  or  equal  to  i , the  instantaneous 

frequency  error  does  not  exceed  the  linear  operating  region  of  the  discriminator. 

Filtering  of  the  discriminator  output , therefore , will  eliminate  the  zero-mean 

fluctuating  component  Af. , and  tracking  behavior  will  not  be  impaired.  For 

A^  greater  than  i , however , the  excursioixs  in  Af.  exceed  the  linear  region 

of  the  discriminator,  and  the  time-averaged  frequency  output  will  be  biased 

so  that  it  lies  somewhat  closer  to  the  interfering  signal  frequency . Thus , a 

signal-to-interference  ratio  which  produces  A.  = i may  be  considered  the 

threshold  for  successful  tracking . In  this  threshold  region , the  approximations 

in  equations  (2.14)  and  (2.15)  do  not  hold.  Based  on  Beghdady's  [l96lj 

exact  results , this  threshold  occurs  for  (S/I).„  = 2 .5  dB . 

ir 

The  results  of  the  preceding  analyses  may  be  summarized  as  follows: 

a . the  tracking  technique  requires  a signal-to-noise  power 
ratio  in  the  IF  bandwidth  of  at  least  0 dB,-: 

b . white  noise  causes  a fluctuating  frequency  error  Af  of  zero 

I& 

BW 

IF 

mean  and  RMS  value  » 

c . the  tracking  technique  requires  a signal-to-interference  power 
ratio  in  the  IF  bandwidth  of  at  least  2.5  dB;  and 
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B. 


d . C«d»3?ent  Interference  et  frequency  f*  causes  an  error  la 

tracking  frequency  denoted  by  Af^,  of  zero  mean  and  RMS 

'5  -^1 

THE  ANTENNA  SYSTEM  AND  THE  INCIDENT  WAVE 
X . Wave  Fields  and  Polarization 


Consider  the  general  incident  wave  described  in  Equation  (1.2)  and 
Fig.  1.1.  Transforming  to  the  x,y,z  ground  coordinate  system,  the  six 
field  components  are 

£ = E [ (“Sin  a - jS  COSO)  (sin  8)  + (cos  a - jR  sin  a)  (-sln^  cos  8)  ] 

E = E [(“Sln  o - jR  cosa)  (-’Cos  0)  + (cos  a - jS  sin  a)  (-sin  |3  sin  8)  ] 

y ^ 

E = E.  [(cos  a - JR  sin  0!)  (->oos^)  ] (2.18) 

Z X 

= Ej^/t)  [(cos  a -jR  sin  a)  (sin  9)  + (sin  a:  + jR  cosa)  (-sin (3  cos  0) 


H = E /rj  [ (cos  Q!  - jR  sin  a)  (-cos  0)  + (sin  a + jR  cos  a)  (-sin  /3  sin 

y 1 • ' 

H = E /tj  [(slnoi  + jR  cos  0!)  (-C08/3) 

Z X 


S) 


If  an  antenna  system  were  constructed  consisting  of  three  orthogonal  electric 
field  sensors  and,  three  orthogonal  magnetic  field  sensors  to  measure  the  field 
components , the  instantaneous  power  flow  in  the  wave  could  be  calculated , 
theoretically , as  the  Poynting  vector , 


P = E H -E  H 
X y z z y 


P = E H - E H 
y z X X z 

P = E H -E  H . 
z X y y X 
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(2.19) 


Measurement  of  the  elevation  angle  of  arrival  is  of  limited  value  because 

of  the  presence  of  earth-ionosphere  refieetions . To  ffiid  the  hh|(le 

of  arrival , only  the  first  two  components , P and  P ♦ nee^  be  determined. 

X y 

however,  the  equations  still  require  aU  six  field  ebmponehls  to  be  meibured. 
Practical  direction-finders  for  ground-based  use  omit  several  of  the  antennas, 
implying  the  assumption  that  the  omitted  field  components  are  smi^l.  ^ 
Typically , a;  system  consisting  of  crossed  loops  and  a vertical  electric 
. antenna  tvill  be  employed , measuring  , and  11^ . Then 


P = -E  H 
X z y 


P 2 E H 
y z X 


(2.20) 


assuming  either  H is  small , or  both  E and  E are  small;  these  conditions 
z X y 

are  realized  with  a purely  vertical  polarization  (in  the  traditional  sense , 

that  is , electric  field  in  the  plane  of  incidence)  for  the  wave,  The  success  of 

this  system  for  direction-finding  on  grovmd-wave  signals  from  vertically 

polarized  transmitters  supports  this  assumption.  When  ionospherically 

reflected  waves  are  received,  however,  the  polarization  may  not  be  purely 

vertical;  the  estimates  of  P and  P from  the  three  field  components  are 
■ X y • 

incorrect , and  the  error  is  referred  to  as  polarization  error.  Working 
backward  from  the  observed  polarization  error  on  a signal  transmitted  with 
vertical  polarization  from  a known  direction  of  arrival , Horner  [l9&?]  found 
that  the  received  polarization  for  a 100-km  path  during  skywave  conditions 


\this0m  and  ea^fionsf . 
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du!^  exit  point  at  lower*  edge  of  the  ionosphaiye , 
idndsphe;^  wav^^ide t to  sitp.  His  ira^tHts  for 

the  polfio^ation  fisp.shG#  wide  ^ a ^<$ibh  pffi^pex^^;  ^ j^ath 

.distance . In  shidyhig  his  plots  * however , the  gen^^  pbsenration  cim  he 
made  that  the  polarization  ellipse  will  tend  to  have  an  axis  ratio  of  about  i , 
and  that  the  handedness  of  the  whistler  polarization  in  the  duct  is  preserved. 
The  conclusion  drawn  is  that  the  direction-finder  design  selected  should  he 
one  which  produces  the  minimum  possible  polarization  error  over  the  entire 
range  of  received  polariaation . 

2 . Antenna  Elements  for  Field  Component  Measurements . 

At  very  low  frequencies , the  antennas  which  can  he  practically  con- 
structed for  direction-Hhding  are  electrically  very  small  compared  With  the 
wavelengths  involved.  While  the  efhciency  of  such  antennas  may  he  low , the 
radiation  patterns  of  small  dipoles  and  loop  are  theoretically  pet'fect . When 
used  as  a receiving  antenna,  the  small  dipole  is  sensitive  only  to  the 
electric  field  along  the  dipole  axis;  the  small  loop  is  sensitive  to  the  magnetic 
field  perpendicular  to  the  plane  of  the  loop . 


A jocmsi^riible  amount  of  cQn^ 

ndative  merits  of  loops  *»n4  dipoles  for  VIJS  reception.  One  consideration  is 

sensitivity  of  the  receiving  system  , which  will  depend  pn  the  efficiency  with 

which  the  antenna  extracts  energy  from  the  incident  >waye  and  couples  it  to 

the  preamplifier.  Papers  clarifying  this  controversy  have  been  published  by 

Wheeler  [l946]  end  Franks  [l964]  . In  theory , the  avaiinble  power  from  any 

small  antenna,  loop  or  dipole , of  any  sise  less  than  is  the  same  for  a 

given  incident  wave  power  density . In  practice , this  aihonnt  of  power  is  never 

extracted  due  to  the  difficulty  of  presenting  the  antenna  with  a cox^ugete-* 

matched  load  impedance  and  due  to  losses  in  the  antenna  and  matching  circuit 

conductors . The  impedance  of  the  small  antenna  is  primarily  reactive , with 

- 1 

the  radiation  resistance  representing  only  a small  part  of  the  impedance.  A 
useful  figure  of  merit  suggested  by  Wheeler  is  the  quantity  p , where 


R 


rad 


ant 


(2,21) 


which  is  the  ratio  of  the  radiation  resistance  to  the  antenna  reactanoe.  It  is 
indicative  of  the  antenna  efficiency  which  can  be  achieve^  with  practical 
coupling  techniques . In  Figs.  2 . 7 and  2 , 8 , the  impedance  and  efficiency 
characteristics  of  small  dipoles  and  loops  as  a function  of  their  size  are 
described.  Most  interestingly . the  value  of  p is  identical  for  both  antennas 
when  the  dipole  length  is  equal  to  the  loc^  diameter. 

Another  consideration  which  has  bemx  advanced  is  the  relative  Immunity 
to  noise  pickup  of  the  two  types  of  antennas . "Noise  pickup"  may  be  of  two 
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types:  first , nearby  sources  of  noise  may  capacltively  or  inductively 
couple  undesired  energy  into  the  antenna;  second,  there  may  be  direct 
collisions  of  charged  particles  with  the  antenna  element . Practical 
experience  indicates  that  the  electric  field  antenna,  or  dipole,  is  more 
severely  ejected  by  both  types  of  noise  pickup . In  particular,  the  author 
has  witnessed  incidents  in  which  a cloud  passing  directly  overhead  caused 
very  high  potentials  to  develop  across  a vertical  electric  antenna’s  terminals , 
causing  arcing  and  rendering  its  output  useless , while  loop  antennas  at  the 
same  location  were  unaffected.  There  are  possibilities  for  controlling  noise , 
both  by  shielding  of  manmade  sources  and  by  providing  weakly  conducting 
dielectric  shields  to  reduce  particle  impingement . 

3,  Antenna  Axis  Perturbations  and  Site  Effects . 

The  sensitive  axis  of  a dipole  or  loop  may  differ  from  the  intended 
direction  for  mechanical  or  electrical  reasons . Mechanical  misalignment  of 
the  antennas  produces  a predictable  cha^e:  the  direction-finding  errors  so 
produced  are  analysed  in  Section  Il.C .4.  Electrical  causes  of  axis 
disturbance  are  more  subtle  in  nature . All  conductors  in  the  vicinity  of  the 
antennas  respond  to  the  incident  waves  also . crewing  electric  and  magnetic 
fields  which  can  couple  into  tn©  antennas  ct^acitively  or  inductively . Ax^ 
imbalance  in  the  construction  of  the  antennas  and  preamplifiers  reduces 
the  Inununity  of  the  ant^ma  to  such  stray  potentials  and  currents . 

Measurement  of  the  incident  electric  field  preoeiUs  particular  problems 
due  to  the  Infiuence  of  the  cables  used  to  carry  power  to  the  antenna 


preamplifiers  or  signals  awe/  from  the  mitennas.  The  dipoles  must  be 
perfectly  balanced  and  the  cables  led  away  orthogonally , or  else  the  potential 
induced  on  the  cable  bundle  by  the  incident  wave  will  couple  into  the  dipole 
causing  a response  to  an  undesired  component  (£  in  Fig.  2.9a)  of  the  field. 
For  the  vertical  electrical  dipole  measuring  , this  balance  is  almost  impos- 
ible  to  achieve  in  practice.  As  illustrated  in  Fig.  2.9(a)  the  capacitance 
between  the  lower  dipole  arm  and  groupd  will  unavoidably  be  greater  than 
that  between  the  upper  arm  and  ground;  as  a result , the  charge  distribution 
induced  on  the  long  feedline  cable  bundle  by  the  E field  component  shewn 

X 

will  induce  a greater  potential  on  the  upper  arm  than  the  lower  arm , and  the 
difference  will  appear  as  a voltage  across  the  dipole  terminals . Another 

approach  is  the  ”monopole”  antenna  design  shown  in  Fig.  2.9(b) . Since 

« 

the  grounding  process  is  imperfect  at  best . the  imaging  plane  displays  a 
much  lower  resistani^  along  the  direction  of  tiie  feedllne  than  in  other  redial 
directions;  as  a result  the  monopole's  image  is  skewed  and  the  equivalent 
axis  of  the  dipole  is  tilted , again  making  it  sensitlye  to  E^.  These  two 
approaches  make  up  most  existing  antenna  installations  for  VI/F  t^servations 
of  electric  fields . Although  the  installation  illustrated  in  Fig.  2.9(c)  would 
be  expected  to  reduce  the  problem,  practical  considerations  dictate  the 
desirability  of  a system  which  can  operate  with  existing  antenna  installations. 
Consideration  of  the  problems  introduced  by  imperfect  measurement  of 
is  given  in  Section  H . C . 5 . 


FIGURE  2.9  VERTICAt  ELECTRIC  ANTENNA  APPROACHES  AND  PROBLEMS . 

(a)  Sieved  ve£<tic^  dipole/  Unbalance  due  to  presence  of 
g^und  causes  coupling  to  fCedUne  potential , and  sensitivity  to  E_ . <b)  GitjUnd 

bCsed  ntonopdev  Finite  ground  resistwce  causes  lowest  resistance  path  to 
exist  along  feedline  shield,  tilting  dipole  equivalent  exis.  causing  sensitivity 
to  E . (c)  Monopole  with  complete  symmetry  maintained  in  grounding  and 

nea^y  objects . Power  lines  will  destroy  symmetry . Sensitive  only  to  E^. 
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One  more  source  of  error  in  the  direction-finder  is  ”site”  error.  These 
errors  are  not  really  part  of  the  instrumentation , as  they  represent 
deviations  in  the  direction  of  propagation  of  the  wave , and  even  a theoretically 
perfect  direction-finder  would  be  affected.  One  of  the  more  familiar  forms 
of  error  is  the  coastal  refraction  effect , which  is  associated  with  the  dis- 
continuity between  highly  conducting  seawater  and  poorly  conducting  ground. 
Kukes  and  Starlk  |^1962  j state  that  this  form  of  error  is  less  than  1^  at 
wavelengths  greater  than  3000  meters  (100  kHs) . Other  possible  sources 
of  error  are  inhomogeneitiea  in  the  earth  conductivity  and  sloping  ground. 
Inasmuch  as  these  errors  are  dependent  on  the  site,  the  most  appropriate 
solution  is  to  perform  a calibration  of  the  site  using  emitters  at  known 
bearing  angles , determining  the  measurement  errors , and  flitting  a curve  to 
the  error  as  a function  of  angle  . 


C.  DIRECTION-FINDING  SIGNAL  PROCESSING 
1.  Analysis  of  New  Technique , 

The  new  technique  for  determining  direction  of  arrival  from  the  signals 
from  the  small-aperture  antenna  system  consists  of  computing  the  time  ■ 
averaged  Poynting  vector  components  in  the  x-  and  y-directicins , m-J  re- 

f 

solving  the  direction  of  arrival  as  the  reverse  of  the  Poynting  /fee  >>;; 

Stated  mathematically . the  Poynting  vector  estimates  are  ' 


^The  reversal  is  due  to  the  fact  the  direction  of  arrival  is  taken  to  mean  the 
bearing  to  the  signal  source  ridher  than  the  direction  of  signal  propagation. 
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(2.22) 


P = -E  H 
X z y 


P = E H 
y z X, 


(where  the  symbol  ” denotes  an  estimated  quantity) , and  the  time-averaged 
values  of  the  Poynting  vectors  are  therefore 


X 


7[E  H*  + E*H  ] P = 7[E  H*  + E*H  1 , (2.23) 

A XT  T XT  ■'  XT  A ^ IT  Xt  IT  V-J  * \ • C 


z y 


z y 


Z X 


Z X' 


From  these  estimates  of  the  average  power  flow  in  the  x-  and  y-directions , 
the  direction  of  arrival  is  the  reverse  of  the  Poynting  vector,  or^ 


(2.24) 


This  new  technique  differs  from  existing  VLI'  direction-finding  methods 
because  nonlinear  operations , that  is , cross-multiplications  between  the 
antenna  signals , are  performed . This  is  made  possible  by  the  real  time 
multichannel  Altering  action  provided  by  the  frequency-tracking  loop; 
without  prefUtering.  such  nonlinear  operations  would  produce  massive 
intemodulation  between  all  the  spectral  components  present . obscuring  the 
desired  information  unless  an  extremely  high  signal-to-noiae  ratio  were 
present . 


^The  use  of  the  imaginary  operator  *’j”  here  is  not  to  be  confused  with  Us 
usual  application  in  the  representation  of  sinusoidally  time-varying  quantities. 
In  equation  (2,24) , the  complex  quantity  -f  - jt*  represents  a vector  In 
space;  the  positive  real  axis  coincides  with  t^e  x ^orth)  direction  and  the 
positive  imaginary  axis  lies  along  the  y (east)  direction.  Complex  notation 
is  a mathmnaticaUy  convenient  means  for  cmnposing  a ve<Uor  from  two 
orthopjnal  components;  the  angle  of  the  vector  tv  ative  to  the  reference  axis 
is  given  by  taking  the  complex  argument  (z — ) , This  representatioi^is 
superior  to  the  conventional  arctangent  function  which  contains  a 180^ 
ambiguity  and  which  can  only  be  used  to  represent  angles  over  a i30  range 
unless  a special  note  is  made  to  examine  the  signs  of  the  numerator  and 
denominator  forming  the  argument  of  the  fonction. 
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The  implementation  of  the  signal  processing  technique  will  now  be 
discussed.  Since  all  three  input  signals  are  mixed  with  a common  local 
oscillator  signal . the  phase  and  amplitude  relationships  of  the  signals  are 
preserved  during  the  translation  to  the  intermedi^e  frequency . The  IF 
voltages  are  denoted  by  V^,  and  V^.  Since  . as  discussed  in  Section 
n.B.2t  the  efficiencies  of  comparably  sized  loops  and  dipoles  are  equal, 
these  voltages  must  be  given  by 


V = E 
1 z 


v„  = m 

2 X 


V =nH  (e  assumed) 
jf  y 


(2.25) 


The  orosa-raxiitipUcations  and  eLigle  resolution  are  performed  in  a novel 

mannel  In  the  hardware . By  applying  the  channel  I (K  ) signal  to  an 

2 

amplitude  Uraiter,  the  phase  information  is  extracted  in  binary  form;  the 
binary  limiter  output  is  then  digitalis  multiplied  with  a binary  "scanning" 
signal  before  application  to  controlling  inversitm  awitches  in  the  two  loop 
signal  channels . A discussion  of  the  hardware  in  detail  is  I'eserved  for 
$e<^n  lil;  for  purposes  of  analysis , the  process  acn  be  modeled  as  shown  in 
Pig.  2.10(a)  and  (b)."^ 

The  cPoas-multlpUcatlon  operations  in  Fig.  2. 10(a)  yield  the  instan- 
taneous Poynting  vector  estimates:  low -pass  filtering  removes  the  twice- IF 
coo^ponent,  Thus, 

fThe  effects  of  the  amplitude  limiter  will  bo  c^sidered,  however,  when  noise 
or  interference  present . 
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FIGURE  2 . 10  MODEL  OF  SIGNAL  PROCESSING  PERFORMED  TO  COMPUTE 

THE  DIRECTION  OF  ARRIVAL  FROM  THE  WAVE  MEASUREMENTS, 
(a)  Model  of  signal  processing  to  compute  the  time- averaged  Poynting  vector 
components  from  the  antenna  signals . (b)  Model  of  signal  processing  to 

resolve  the  direction  of  arrival  angle  from  the  time- averaged  Poynting 
sector  components . 
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(2.26: 


‘[V2  *^'^2]' 


These  voltages  are  then  multiplied  by  the  scanning  signals  as  in  Fig.  2.10(b) 
and  summed  to  produce  a signal  of 


A jW  t ..  jW  t 

'IT  - ^ ®can  — scan 

Ve  - -»)P^e  -jT?P^e 


(2.27) 


This  voltage  is  compared  in  phase  to  the  scanning  carrier  e 


jw  t 


Then, 


A / A J<a  t A t / 

e . /iP^e  -jnPye  “““  - 


Which  is  equivalent  to 


8 = /-T)P  -jt)P  + W t - t 

/ X y scaa  scao 


thuS(  the  desired  relation  of  (2.24)  has  been  implemented , 


e = Ap  -jp 

/ X •’  y 


To  demonstrate  the  validity  of  the  computsiitms , let  us  temporarily 
assume  an  incident  wave  at  angle  6 , purely  vertically  polarized  so  that  a = 0 
and  R = 0.  The  field  quantities  in  equations  (2.18)  now  reduce  to 


cosp 


Ej/n  sin© 


(2.29 


and  the  IF  voltages  are 


-Ej^/n  cos© 


-Ej^  cos  p 


Ej  sin© 


-Ej  cos© . 


■ : :;s: 


The  cross-multipUgetions  produce 


cos ^0080 

(2.31) 

— 2 
’IPy  «=  "i  cospsine 


end  the  angle  resolution  produces 


or 


9 = 


cosp  (CO80+  j sin0  ), 


(2.32) 


0=0 


(2.38) 


The  technique  provides  a unique , unambiguous  indication  of  the  direction 
of  arrival . The  causes  of  error  in  5 will  now  be  investigated . 

2.  Errors  Due  to  Noise  and  Interference . 

The  ejects  of  internally  generated  noise  in  the  equipment  provide  an 
ultimate  Umiiation  to  the  direction-finder  sensitivity  for  weak  signals . The 
total  noise  contributions  in  each  channel  may  be  lumped  into  an  additive  noise 
voltage , denoted  by  n,  (t) , n*.  (t) , and  n, (t)  for  the  three  channels . These 
three  voltages  are  Independent . originating  in  separate  processes , but  they 
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may  js^pe  ccmside««td  to  have  Gaussian  statistics , with  sero  mean,  and 
idend<B^  mean-squa^d  values: 


«.(t)  = n (t) 

A 


ng(t)  = 0 


(2.34) 


2 2 2 2 
nj^"(t).  = n^(x)  = (t)  = n^(t) 


To  analyze  the  effects  of  these  noise  voltages , we  can  adopt  a.  model  for  narrow 
band  noise  presented  by  Panter  j^l965j  containing  independent  Gaussian 
components  in  quadrature  at  the  center  frequency  of  the  band.  The  noise 
voltages  at  the  IF  circuitry  are  thus  represented  by 


n(t) 


n (t)e  + jn^  (t)e 
C 8 


(2.35) 


where  n^(t)  and  n^(t)  are  slowly  varying  compared  to  e 


and  for  which 


n^(t) 

c 


n^(t)  . = n^t). 

s 


(2.36) 


It  is  desirable  to  express  n (t)  in  terms  of  the  signal-to-noise  power  ratio  in 
the  IF . To  avoid  ambiguity  due  to  the  fact  that  the  signal  varies  with  the 
incident  wave  geometry , we  will  define  (S/N)^^  to  be  the  signal-to-noise  ratio 
which  would  exist  on  a comparable , but  nondirectional , receiving  system . 
Then. 


V 


iw  t 
•'if 

(e  assumed) 


and  the  signal  power  is 
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The  ©oiae  ptw^is  gtyen  hy;  n ; (t) , laia^^ere&isa  v; , v ; v-  ; 
(s/N)jj.  * z^nnity  r 


n^(t)  = Ej^/2(S/K)jp  (2.37) 


(Hereafter,  the  time-dependence  of  n , n » n will  be  implicit.)  The  IF 

c s 


(e^^F^  assumed  for  all  quantities) 


0 

and , for  vertical  polarization  of  the  incident  wave , ( a = 0 , R - 0) , we  have 

= -Ej^cosp  + n^, 

= E^sin9  + n^2  (2.39) 

= -Ej  0088  *1.^3  * in^3 


At  this  point  we  must  taken  into  ap^ont  the  that  a limiter  amplifier  is 
included  in  channel  1.  As  discussed  in  Section  II. A. 2,  the  noise  is  reduced  in 
the  presence  of  the  signal  when  (S/N)^  > 1 , aocordingto  the  of  (2 . 10) , 


(S/N) 


LIM 


2 (S/N) 


IF" 


We  must  therefore  modify  n^j^  and  n^^  in  by  the  factor  l/'s/2.  Proceeding 

with  the  cross-multiplications  of  Vj^  with  and  with  , and  the  low -pass 

filtering,  we  derive 


r?P 


X 


1 2 

-r  [E,  008  fl  008  6 -(1/ nT2)  E.  cose  -E,  cos/ja.„ 
2 1 ' i ui  i.  oo 


.)P  = I [-Ej^  008  3 sin  8 +(l/'/2)  Ej  sln0  -E^  ooBJn^j 


(2.40) 


Besolving  the  direction  of  arrival  using  equation  (2.24)  we  have,  after  some 
rearranging, 


(Ej^  008  0 - ^ E^  (008  e + ) sUx  g) 


* < 4 "ol  ““  « %3  -J‘o2>  *d-  M <“83  -J“82» 


(2.41) 


» i • 


Expressing  this  as 


Be*®  * 


and  dividing  by  Be^  ^ (assuming  B > 0) 


(A/B)e^  = I + (C/B)e^  , 


A 

the  angle  error  A6s  0-0  is  given  by 


A0-  /l+(C/B)e^’'^^  = arctm  (C/B) 


arctan 


n , n „ -n  . n „ 
ul  u7  u4  tt9 


-2ii  -J2  n , 
g6  ^ ul 


oosfl  (8/Nl 

”u4  \s 
3^/2  coe^  (S/ISf) 


From  this  expression  we  can  see  that  the  error  has  a zero  mean.  To  deter- 
mine the  RMS  error,  a <x>mputer  slmulati<xi  v;as  performed  using  randota 
nutcber  generators  shaped  to  have  Gaussian  distidbutions.  These  reetUts  are 
shown  in  Figure  ZAl.  For  hifi^  signal-tO'*nolse  ratio  oonditions , the  error 
can  be  approximated  by: 
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FIC€RE  2 . 11  RMS  DIRECTION-FINDING  ERROR  VERSUS  SIGNAL-TO-NOISE 

RATIO  WHICH  WOULD  BE  PRODUCED  ON  A COMPARABLE , BU 
NONDIRECTIONAL,  RECEIVING  SYSTEM- 


and  therefore 


iTOTN) 


We  will  now  consider  instead  the  effects  of  a second  coherent  signal 
within  the  receiver  passband.  The  ratio  of  the  desired  signal  power  to  the 
interfering  signal  power  is  d'eiioted  by  (8/1)^^.,  and  the  difference  frequency 
between  the  two  signals  by  f^ . Whereas  the  azimuthal  angle  of  arrival  of  the 
desired  signal  is  0 , the  angle  of  arrivsl'.of  the  interference  wlU  be  6 
the  elevation  angle  of  the  desired  signal  is  p t and  the  elevaffon  angle  of  the 
interfering  signal  is  (3^ . 

The  asipUtude  of  the  interfering  signal  oan  be  oeloulated  as: 


E,/  v?s75;7 


and  the  signals  in  the  IF  are  therefore  given  by  (e  assumed  for  all 
quantities) 


“ 1W5Z  ° 

IJ? 

sin  (0+&)  jw  t 

(2.60) 

lx 

cos  (6+^)  jw  .t 

= -Et  °°B  6 -E^  :jmr'  ° 

lx 

At  this  point  we  must  again  take  into  account  the  presence  of  the  limiter 
ampliCier  included  in  the  cross-multiplication  circuitry . Incorporating  the 
S/N  enhancement  from  equation  (2. IS)  and  performing  the  cross^ultiplications , 
we  have 


Tjp  » - 


cos  w ,t 

2 ^2  u 

cos  ^ 008  8 + Ej  008  ^ 008  (0+0^)  • • 

IF 


0 2 
+ E^“  cos  0^  cos  0 ‘ * ’ 

IF 


2 (S/I) 


IF 


(2.51) 
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'-!5 

??  ’u 


TIP  ^ 

y 


1 r 2 2 

~2  \'h  COB  P slniB^B^) 


2 '-^d*  2 

..  -E^  \ cos^^sin(0+e^) 

IF 


••*2  (S/I) 


1 1 

A)_  * 


IF  J 


The  resolution  of  these  quantities  into  an  estimated  angle  6 gives 


Exptt'ssirg  this  as 


J0 


cosp 


and  dividing  by  Be''  . (assuming  < 1) 


(A/B)e 


j(9-8)  ^ 


(C/B)e 


J8d  . 


thus , the  angle  error  A6  = e - 8 , is  given  by 


A0 


= jx^ 


(C/B)e 


iGd 


(2.54) 


(2.55) 


or 


(2.56) 


From  this  expression  we  can  see  that  the  error  contains  both  a non-zero  mean 
value  and  a periodic  variation  at  the  difference  frequency  . The  character 
of  A0  versus  time  is  Illustrated  in  Fig.  3.12  for  some  representative 
conditions . The  mean  error  A8  and  the  RMS  varieiion  about  the  mean  for 
various  cionditlotm  are  plotted  in  Fig.  2.13. 

Approxlmationa  can  edso  be  made  \ ^ich  provide  a simpler  indication 
of  the  effects  of  interference.  For  ^h  signal -lo-inierference  ratio 
conditions*  the  error  is  well  approximated  by 
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FIGURE  2.13  MEAN  AND  RMS  DIRECTION-FINDING  ERROR  VERSUS  (a)  SIGNAL- 
TO-INTERFERENCE  RATIO;  Cb)  ELEVATION  ANGLE  RATIO: 

(c)  RELATIVE  INTERFERING  SIGNAL  BEARING  . 


A6 


cos  ^ 


cos  Wjt 
cl 


2 cos  (i  (S/I)  •/(S/l) 


sine. 


(2.57) 


demonstrating  that , 


_ cos 

""  2 cos  /3  (S/I) 


(2. 58) 


and 


^®BMS  ~ sf2^/l)  ' ®d 


(2.  59) 


3,  Polarization  Error 

To  this  pcdnt  in  the  analysis  we  have  restricted  the  polarization 
of  the  incident  wavs  to  purely  vertical . When  this  condition  is  relaxed , the 

A 

estimates  of  the  Poynttng  vector  components  and  P^  based  only  on 
measurement  of  t , and  may  be  incorrect , and  the  computed  angle  of 
arrival  may  contain  an  error  known  as  polarization  error. 

Let  a and  R take  on  nonzero  values  so  that  a general  polarization  Is 
permitted.  The  fields  from  equations  (2.18)  are  measured  by  the  antennas, 

jwjjft 

and  presented  as  IF  voltages  (e  assumed), 
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(cos  a - jR  sin  o!)  (-cos  j3) 

^2  ” Q!  ~ jR  sin  a)  (sin  6)  + (sin  a + jR  cos  a)  (-sin  cos  B) 

Vg  = (cos  0!  - jli  sin  a)  (-cos  0)  + (sin  a + JR  cos  a)  (-sin^  sla0). 

The  cross-multiplications  then  yield 

A 

~ 12  r 2 2 2 

= 'i  a+R  sin  ) (cos  0)  . . . 

2 1 
. . . + (1-R  ) (cos  a sin  a)  (sin  /3  sin  0)J 

(2. 

1 2 r 2 2 2 

^Py  = 2 COS ^ [(cos  0!  + R sin  a)  (-sin  0)  . . 

2 -1 
...  + (1-R  ) (COSO!  sin  O')  (sin^  cos  0)1 

Resolving  this  into  an  angle  8 , 

/ 


e 


[(cos^  0!  + R^  sin^  + [sin^  (R^  -1)  cos  a sin  a] 


and  expressing  this  as 


a/«  = Bei®  . 

l6 

we  divide  by  Be’'  (B  ^ 0 , providing  a polarization  other  than  pure  linear 
horizontal  is  assumed) , 

(A/B)  ^ 1 + (C/B)J'^^^ 


or 


Ad 


y'i  + (C/B)  = 


arctan  (C/B) 


so  the  polarization  error  is  given  by 


A6 


arctan 


r 2 

sin  B (B  -1)  cos  a sin  a 

2 2 2 
L cos  a + B sin  a 


(2.62) 


The  polarization  error  as  a function  of  R,  p,  and  a is  plotted  in  Fig.  2.14. 


For  purposes  of  comparison , the  polarization  error  of  the  goniometer  system, 


based  on  the  analysis  of  Bullcugh  and  Sagredo  1973  , is  plotted  in 


adjacent  frames.  The  advantage  of  the  new  technique  is  evidenced  by  lower 
errors,  especially  for  high  elevation  angles  and  polarizations  near  circular. 


aCdegs) 


FIGURE  2 . 14(a) . POLARIZATION  ERROR  FOR  NEW  TECHNIQUE  AND 

conventional  goniometer  as  a function  of 

ELEVATION  ANGLE  AND  POLARIZATION  ELLIPSE  AXIS  TILT  FOR  R « 0 


a (deqs) 


FIGURE  2 .14(b) . POLARIZATION  ERROR  FOR  NEW  TECHNIQUE  AND 
CONVENTIONAL  GONIOMETER  AS  A FUNCTION  OF 
ELEVATION  ANGLE  AND  POLARIZATION  ELLIPSE  AXIS  TILT,  FOR  R i 
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FIGURE  2 . 14  (d) . POLARIZATION  ERROR  FOR  NEW  TECIiNlQUE  AND 
CONVEtmONAL  GONIOMETER  AS  A FUNCTION  OF 
ELEVATION  ANGLE  AND  POLARIZATION  ELUPSE  AXIS  TILT . FOR  R = 3/4 
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Antenna  Misalignment  Errors . 

In  Section  H.B.3  it  was  pointed  out  that  the  sensitive  axes  of 

the  antennas  may  differ  from  the  desired  directions . Consider  a system  of 

antennas  where  the  antennas  measuring  H and  H are  misaligned , because 

X y 

of  me<dianioal  or  electrical  asymmetries » so  that  the  sensitive  axes  are 
perturt>ed.  Referring  to  Fig.  2 . 15 » the  projections  of  the  misaligned  axes  on 
the  xy  plane  are  displaced  by  angles  5 and  ( (in  radians) . Now  the  fields 
measured  by  the  antennas  are  in  error,  as  follows: 

A 

H = H cos  6 + H sin  6 
X X y 

(2.68) 

H = H eoB€  - H Sine  . 

y y X 

Substituting  the  actual  fields  from  equations  (2.18)»  and  assuming  ns  0 and 
E 5=  0 for  simplicity , the  IF  voltages  will  be  given  by 

cos  p 

V =*  E sin  9 cos  5 -E  cos  6 sin  6 

4W  .L  A 

(2.64) 

V = *E  cos  6co0  c ~E  Sin  9sln  f . 

•J  i.  ■ • 

The  cross- mMUiplicatlons  then  yield, 
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A 


2 

“iEj^  (cosjJ  cose  cos€  + Qosp  sine'  sine) 

2 

iEj^  (“cosp  sine  cos  8 + cosp  cos0  sin6) 


Resolving  the  angle,  and  assuming  8 and  c are  small,  cos 6 , cos  t 


6 = 


Denote  two  new  angles  and  x as  follows: 


+ sint  sine  jsin6  cose 


= 1 , 

(2.65) 


^ = ( « + 6)  /2 
X = (‘  -8)  /2 


(2.66) 


Note  that  4»  is  the  rotation  of  the  axis  of  symmetry  of  the  antennas  due  to  the 
misalignment,  and  X is  the  misalignment  of  the  y-axis  antenna  relative  to  the 
new  axis  of  symmetry . With  these  quantities , 


Q - 


+ sin  (x+t/;)  sin  6 + j sin  (x~0)  cos  $ 


and  using  trigonometric  identities . 


e 


+ slnx  COB  (|)  ]e  - sin  ^ cosX  je^® 
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Since  4<  and  X are  small 


B - “ 5)  ^ .M-./Z) 


$ & 


(2.67) 


Expressing  this  a*i 


je 


eJ6f  = e-’'^  + jXe'^"  - jve' 


jB 


and  dividing  b 


1 + -j.4^  . 


we  can  see  that  for  X and  4^  small , 

A 

B-S  ^ X cos  2S  - 4* 
or 


A0  = - ^ + X 


(2.68) 


The  mean  value  corresponds  to  the  cdsaligttment  of  the  axis  of  symmet  ry . 

The  periodic  (qui^ran(ai)  c'j.tipotieni  has  a peak  value  corresponding  to  the 

misalignment  of  the  y-ajds  isnteinia  relative  to  the  axis  of  symmetry* 

Now  consider  a ndsaiigriiaent  of  the  loop  antennas  corresponding  to 

non  hotiz<mtality  of  the  equivalent  axes  of  the  loops , as  shown  in  Fig  2.16. 

> 
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y 

(EAST) 


FIGURE  2.16.  ISOMETRIC  VIEW  OF  DIRECTION  FINDER  SITE  SHOWING 
TILT  IN  SENSITIVE  AXES  OF  LOOP  ANTENNAS . 


Let  the  tilt  angles  of  the  and  loops  be  denoted  by  and  (in 
radians) . Then  the  Poynting  vector  estimates  will  again  be  in  error , as 

Py  = \ [“,““’•2  - 

Assuming  the  tilts  are  small , and  using  the  field  quantities  from  equation 
(2.29)  for  0'  = R = 0,  the  resolved  angle  is  given  by 

A 

e = 

or , equivalently , 


AS  - Tg  cos  8 - Tg  sin  0 . 


(2.69) 


Finally , consider  the  alignment  of  the  vertical  antenna.  In  the  new 
direction-finding  signal  processing  technique  presented  herein , a small 
error  in  the  measurement  of  does  not  produce  any  first-order  effect  on  the 
direction-finding  error , since  the  mcaiurement  discrepancy  simply  appears 

A A 


as  a common  factor  in  P and  . The  polarization  error  characteristics  will 

X y 

be  slightly  modified  but  only  in  nature  and  not  overall  magnitude . This 
is  fortunate , as  it  was  pointed  out  in  Section  II.B .3  that  the  vertical  electric 
dipole  is  especially  susceptible  to  axis  perturbation.  However, 
tecliniques  have  been  proposed  which  depend  critically  on  the  vertical 
Vienna  performance . and  this  is  the  topic  of  the  next  section . 
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5. 


Errors  in  "Quadrature  Direction-Finders”  Due  to  Antenna  Axis 


Perturbation 

Cousins  1^1972  and  Tsuruda  and  Hayashi  |^1974  j among  others , 
have  proposed  and/or  implemented  VLF  direction-finder  systems  designed 
to  have  zero  polarization  error.  They  are  based  on  utilization  of  the  loop 
signal  components  which  are  in  exact  quadrature  with  the  signal  on  the 
vertical  antenna , as  opposed  to  the  present  system  which  uses  the  in-phase 
component . Recalling  the  eiq)re88ions  for  the  IF  voltages , they  can  be  recast 
as 


[a  ( -cos  (3)] 

V - E !’A(sine)  + B (-sin  |3  cos  6)  ] (2.70) 

2 ^ 

V = E [a  (-cos  6)  + B (-sla/3  aiti  0)  ] 

3 i 

where; 


A = cos  a - JRsId  a 
B = sin  a + jRcos  o 


(2.71) 


The  B-coroponents  are  not  wanted  in  the  in-phase  technique.  That  portion  of 
B which  is  in  phase  with  A,  that  is,  A • B,  causes  polarization  error.  By 
contrast , in  the  quadrature  technique , the  A components  are  rejected  by 
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shifting  the  signal  phase  by  90°  before  the  cross -multiplications . With 
only  the  B components  left , the  angle  is  resolved  in  a different  way . 
Representing  the  operations  mathematically , 


I FJVj  V*  + Vj  ] + I [jv^  v’  + oVj)*  V3  ] 


(2.72) 


or 


e'  = 


/r  • ’ 

/JI  AB  -A  B]  sin  j3  cos  )3  (cos  6 + jsin  0) 


(2.73) 


yielding 


e' 


E 


1 lA 

R cos  jS  sin  ^ e 


so 


e'  = 9 . 


The  quadrature  direction  finder  displays  no  polarisation  error  according  to 
the  formula  shown . However , note  that  the  quadrature  direction  finder 
does  not  produce  a dii'ectional  indication  for  linear  polarization  (R  = 0) , or 
if  0°  or  90°;  whereas  the  in-phase  technique  only  fails  if  the  polarization 
is  pure  .‘inear  horizontal  (R  « 0 and  a-  s 90°  simultaneously) , or  if  |3=  90°. 
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Now  introduce  a aUght  tilt  in  the  sensitive  axis  of  the  vertical  electric 
antenna , by  an  angle  t (in  radians) . The  tilt  will  be  assumed  to  be  toward 
the  x-axis , for  simplicity , but  the  analysis  is  easily  extended  for  general 
titls . The  erroneous  fijrid  measurement  is 

= E + rE  , 
z z x 

The  measured  voltage  is  now 

Vf  = [a  (-COS/3)  + rB  (-sin  6)  + t A (-sin  jj  cos  B)  ] 

The  resolved  angle  is  now  modified  by  the  addition  of  the  extraneous 
component , so  that 


IE 


“ /y  [b  cos/3  sin/5  e-*®  + tR  sin^ /3  cos 


...  + rR  (cos^  a - sln^  a)  sin  8 


(2.74) 


The  error  under  this  condition , for  t small , is  given  by 


AS'  = -r  [■ 


2 cos  (2a)  sin  9 
sin  2^ 


(2.76) 
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The  directional  error  is  slnusoidd  with  aalmuth , with  the  maximum 
error  occurring  at  assimuths  perpendicular  to  the  direction  toward  which  the 
antenna  is  tilted.  The  peak  directional  error  in  degrees  per  degree  of  axis 
tilt  is  given  by  the  ratio  A97Tsin6,  which  is  picked  in  Fig.  2.17.  The 
sensitivity  of  the  system  to  vertical  antenna  orientation  and  perturbation  is 
demonstrated  by  the  high  valuea  of  this  ratio. 
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TILT  ANGLE  OF  POLARIZATION  ELLIPSE  a 


FIGURE  2.17.  DIRECTION  OF  ARRIVAL  ERROR  IN  THE  QUADRATURE 
DIRECTION  FINDER  DUE  TO  A TILT  IN  THE  EQUIVALENT 
AXIS  OF  THE  VERTICAL  ELECTRIC  ANTENNA,  FOR  VARIOUS  VALUES  OF 
ELEVATION  ANGLE  AND  POLARIZATION  OF  THE  INCIDENT  WAVE. 
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m.  DESCRIPTION  OF  THE  APPARATUS 

To  demonstrate  the  feasibility  of  the  new.  ^j>roach  to  VLF  direction- 
finding and  data  gathering  described  in  Fig.  1.3,  a prototype  instrument 
embodying  the  new  principles  and  techniques  presented  was  designed  and 
constructed  by  the  author  at  Stanford  Electronics  Laboratories  . In  this 
chapter,  a description  of  the  apparatus  at  the  block  (hagram  level,  with  some 
insight  into  the  more  novel  circuits , will  be  provided . Detailed  schematic 
diagrams  are  presented  in  Appendix  B. 

A.  SYSTEM  DESIGN 

A diagram  of  the  instrument  emphasizing  the  fun<'tional  blocks  is 
presented  in  Fig.  3.1. 

The  output  voltages  from  the  antenna  block  contain  the  desired  signal 
in  addition  to  the  undesired  energy  in  the  remainder  of  the  VLF  spectrum . 

The  translation  and  filtering  block  serves  first  to  translate  the  input  spectrum 
so  that  the  desired  signal  lies  at  the  intermediate  frequency  (IF);  then  by 
passing  the  translated  spectrum  through  the  IF  bandpass  filter,  a filtered 
version  of  the  desired  signal  (V^e  ) is  obtained.  The  filtered  signals 
from  the  three  antenna  elements  are  applied  to  the  direction-finding  processor 
in  which  the  estimated  direction  of  arrival  is  resolved . 

Turning  to  the  lower  part  of  the  diagram , the  parameter  extraction 
block  makes  measurements  of  the  frequency  error  (Af) , signal  strength  (SS) . 
and  signal -to-noise  ratio  (SNR)  existing  in  the  IF.  Based  on  the  latter  two 
quantities , the  signal  detection  block  makes  a decision  as  to  whether  a 
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FIGURE  3.1.  BLOCK  DLA.GRAM  OF  THE  APPARATUS  . 


desired  signal  is  present  • If  an  affirmative  decision  is  made , the 

firequency  tracking  block  is  enabled.  Under  this  condition , the  frequency 
error  information  from  the  parameter  extraction  block  is  permitted  to  control 
the  tuning  frequency  of  the  receiver  to  accomplish  frequency  tracking  of  the 
signal . 

In  the  data  display  and  recording  block , the  directional  data  and 
additional  parameters  of  interest  are  made  available  for  displ£^  and  chart 
recording . In  addition , the  analog  voltages  are  converted  into  a frequency- 
modulation  format  and  frequency-multiplexed  into  a single  composite  signal 
which  can  be  recorded  on  an  ordinary  audio  tape.  This  approach  facilitates 
the  storage  and  correlation  of  the  data  with  the  raw  VLF  spearum  from  which 
it  was  extracted. 

B.  ANTENNA  SYSTEM 

The  instrument  was  designed  to  operate  with  the  existing  antenna  and 
preamplifier  systems  at  the  Stanford,  California,  and  Roberval,  Quebec, 
observing  sites . The  physical  design  of  the  antenna  system  is  shown  in 
Fig  .3.2.  The  preamplifiers  are  desigiaed  so  that  the  voltage  outputs  from 
the  loops  and  vertical  antennas  are  equal  and  in-phase  for  a given  incident 
signal  (allowing  for  the  directional  patterns  of  the  elements) . These  outputs 
therefore  correspond  to  E , , md  hH  . 
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C.  TRANSLATION  AND  FILTERING 

The  operations  accomplished  within  the  translation  and  filtering  block 
are  illustrated  in  Fig.  3.3.  The  local  oscillator  signal  is  split  and  heterodyned 
against  each  of  the  preamplified  antenna  signals , generating  sum  and  differ- 
ence frequency  components . When  the  local  oscillator  frequency  is  set 
higher  than  the  signal  frequency  by  an  amount  exactly  equal  to  the  inter- 
mediate frequency , the  difference  component  falls  in  the  passbend  of  the  IF 
filters  and  the  signal  is  passed. 

The  bandpass  filtering  is  accomplished  by  active  filter  circuits 
Tow , 1969  j constructed  from  operational  amplifiers . To  provide  a response 
with  a flat  top  and  steep  skirts , two  resonators  are  stagger-tuned  to 
frequencies  just  above  the  below  the  IF  center  frequency.  The  response 
characteristic  obtained  is  shown  in  Fig.  3.4.  The  selection  of  the  IF  filter 
bandwidth  was  the  result  of  a compromise  between  rejection  of  interference 
and  tracking  capability  for  signals  which  change  frequency  rapidly . 

According  to  the  relationship  derived  in  equation  (2.7  ) , the  340  Hz  band- 
width selected  permits  tracking  of  frequency  changes  up  to  57.8  kHz  per 
second . 

D.  DIRECTION-FINDING  PROCESSING 

The  direction-finding  teelinique  has  been  analyzed  in  Section  II.C  .1. 

The  illustration  of  Fig.  2.10  provides  a model  for  the  processing  operation, 
but  in  practice  the  use  of  analog  multipliers  is  avoided  by  applying  the  novel 
implementation  shown  in  Fig.  3.5.  The  IF  signal  from  the  vertical  antenna 
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FIGURE  3.3  OPERATIONS  WITHIN  THE  TRANSLATION  AND  FILTERING 
BLOCK. 
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FIGURE  3.5  . DIRECT?ON-FINDER  PROCESSING  BLOCK  DIAGRAM 


J 


is  applied  to  an  amplitude  limiter  to  convert  it  into  a binary  waveform 
retaining  only  the  phase  information  from  the  signal . Similarly , the  multi- 
phase oscillator  produces  binary  waveforms  corresponding  in  phase  to  the 
jco  t ju  t 

sinusoids  e and  -je  . The  multiplications  can  then  be  perform 

between  the  binary  waveforms  by  using  digital  exclusive-or  gates , and 
multiplication  of  the  binary  waveform  with  the  analog  signals  and 
requires  only  an  inverting/ noninverting  switch . Upon  ba  '.dp, ass  Sitering 

of  the  combined  outp  it  a.t  the  frequency  w , the  higher  harmonics 
generated  by  the  use  of  binary  waveforms  are  removed  and  the  summed 
voltage  corresponds  to 


'"a  “ " ''i  V ^ 


jW  t 
scan 


jw  t 

When  the  phase  of  this  voltage  is  measured  relative  to  e , the  angle 

estimate  is  orodueed , 


fi,  PAStAMETES  EXTRACTION 

The  operational  functions  of  the  parameter  extraction  block  are  shown 
in  Fig.  3.6.  The  filtered  IF  signal  from  the  vertical  antenna  channel  is 
applied  to  a frequency  discriminator  to  measure  the  frequency  error  Af . The 
signal  strength  SS  is  determined  by  applying  the  same  signal  to  a logarithmic 
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FIGURE  3.6  OPERATIONS  WITfflN  THE  PARAMETER  EXTRACTION 


amplifier,  detector,  and  low-pass  filter.  The  use  of  a logarithmic  amplifier 
permits  the  detector  to  operate  on  signals  over  a 60  dB  dynamic  range . 

The  unfiltered  VLF  spectrum  from  the  vertical  antenna  is  applied  to  a 
receiving  subsystem  identical  to  the  local  oscillator , mixer , and  IF  filter 
used  for  the  three  primary  channels . This  receiving  subsystem  is 
independently  tuned , manually , to  a region  of  the  VLF  spectrum  where  the 
signals  of  interest  are  absent;  broadband  noise  originating  in  spheric 
activity , however,  ia  present  in  this  region.  Thus , when  this  independent 
channel  is  applied  to  a logarithmic  amplifier,  detector,  and  low -pass  filter, 
a measure  of  the  noise  amplitude  is  obtained.  By  subtracting  this  logarithmic 
value  from  the  logarithmic  signal  strength , the  logarithmic  signal-to-noise 
ratio  is  calculated . 

F.  SIGNAL  PRESENCE  DETECTION 

The  function  of  this  block  is  to  determine  whether  a signal  of  interest 
is  present . Thresholds  are  set  for  signal  strength  SS  and  signal-to-noise 
ratio  SNR  which  must  be  exceeded  to  cause  the  transition  from  a no-signal 
decision  to  a signal-present  decision . The  reverse  process  can  be  initiated 
by  a loss  of  signal  strength  or  signal-to-noise  ratio , in  addition  to  auxiliary 
checks  on  the  signal  frequency  and  the  time  duration  of  the  signal . 

G.  FREQUENCY  TRACKING 

The  frequency  tracking  block  consists  of  an  integrator  with  a pre- 
settable initial  condition  and  with  electronically  switched  integrate/hold  modes . 


The  input  to  the  integrator  is  the  frequency  error  information  from  the 
parameter  extraction  block . The  output  is  a voltage  which  controls  the 
frequency  of  the  receiver  local  oscillator. 

In  operation , the  operator  selects  the  initial  condition  as  a frequency 
which  the  signal  of  interest  is  likely  to  pass  through . When  the  signal  first 
appears  in  the  receiver  passband  and  is  detected  by  the  signal  presence 
detection  block,  the  integrator  is  enabled.  Thereafter,  any  offset  between 
the  signal  frequency  and  the  center  frequency  of  the  receiver  passband  is 
detected  by  the  parameter  extraction  block  and  applied  to  the  integrator  to 
correct  the  tuning  frequency  of  the  receiver.  When  the  signal  has  terminated, 
the  integrator  is  reset  to  the  initial  condition  frequency  to  await  another 
signal. 

H.  DATA  DISPLAY  AND  RECORDING 

The  direction-of- arrival  data  and  signal  parameters  of  interest 
(frequency . signal  strength , and  signal-to-noise  ratio)  are  made  available 
for  display  on  a panel  meter,  digital  readout,  or  external  analog  chart 
recorder.  In  addition , the  frequency , direction-of- arrival , and  either 
signal  strength  or  signal-to-noise  ratio  are  converted  to  a frequency- 
modulation  format  and  combined  into  a single  audio  channel  suitable  for 
recording  on  a track  adjacent  to  tlie  recording  of  the  broadband  VLF  spectrum 
which  produced  the  data. 

The  tape  data  interface  is  diagrammed  in  Fig.  3.7.  The  local 
oscillator  is  heterodyned  against  an  oscillator  operating  at  the  IF  frequency 
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nOUBE  3.7  CONVERSION  OF  DATA  TO  A FREQUENCY-MODULATION 


to  produce  a difference  ft'equency  component  at  the  signal  ^?equency . The 
sum  hi'equency  component  is  removed  by  a low -pass  filter.  Thus,  the 
tracking  firequency  of  the  receiver  is  represented  by  a sine-wave  at  the 
tracking  frequency  itself.  The  directional  and  signal  strength  data  are 
applied  to  voltage-controlled  osdllators  to  generate  frequency  modulated 
sine  waves  representing  6 and  log  SS . The  directional  data  reside  in  the 
10.4  to  17.6  kHz  band  and  the  signal  strength  data  are  in  the  18.0  to  19.2  kHz 
region , Both  of  these  signals  can  be  combined  with  the  frequency-indicating 
signal , which  is  confined  to  the  0 to  10  kHz  band , by  a simple  addition.  The 
resultant  signal  carries  all  the  information  in  a 0 to  20  kHz  bandwidth . 

In  Fig.  3.8,  an  example  is  presented  showing  a spectrum  analysis  of 
the  data  channel  adjacent  to  a spe*:^ rum  analysis  of  the  raw  VLF  spectrum 
from  which  the  data  were  extracted.  A useful  method  for  applying  the  dsla 
oltannel  is  also  shown , in  which  the  raw  VLF  spectrum  is  mixed  with  the  data 
channel  before  spectrum  analysis . The  tracking  of  the  whistler  is  evidenced 
by  the  superposition  of  the  tracking  frequency  trace  and  the  whistler  trace . 
For  detailed  examination  of  the  direction- finding  and  signal  strength  data, 
a separate  expanded  analysis  of  the  10  to  20  region  of  the  data  track  can 
be  performed  as  will  be  demonstrated  on  the  field  data  records  in  Section  IV. 
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FIGURE  3,8  UTILIZATION  OP  FM  DATA  CHANNEL,  (a)  Conventional 
VLF  spetSiiira  anaiy sis  of  a whistler  received  at  Stanford . 

(b)  Spectrum  analysis  of  the  FM  data  channel  produced  during  reception  of 
the  whistler , Traces  indicate . top  to  bottom , the  relative  signal  strength , 
direction  of  arrival,  and  frequency  of  the  signal,  (c)  Analysis  of  composite 
audio  formed  by  mixing  the  raw  VLF  data  with  the  FM  data  track . Tracking 
of  the  signal  is  evidenced  by  the  auperpositicm  of  the  tracking  frequency  trac« 
and  the  whistler  trace . (d)  Expanded  version  of  (c).  This  record  also 

shows  the  100  Hz  tone  used  to  mark  the  event  for  rapid  access  at  a later  date , 


IV.  EXPERIMENTAL  OBSERVATIONS  AND  RESULTS 


After  initial  testing  of  the  prototype  frequency-tracking  direction- 
finder in  the  laboratory , the  equipment  was  installed  first  at  Stanford , 
California,  and  then  at  the  Stanford-operated  observing  site  at  Rbberval , 
Quebec,  Canada.  These  testing  programs  are  discussed  herein  according 
to  chronological  order.  The  quality  and  degree  of  interest  in  the  results 
actually  follows  the  reverse  order;  partly  because  the  development  of  the 
equipment  was  still  in  progress  during  the  Stanford  program , and  partly 
because  of  the  greater  abundance  of  VLF  signals  of  interest  at  Roberval . 

A.  LABORATORY  TESTING 

1‘  Frequency  Tracking 

After  testing  of  the  individual  circuits  and  fimctional  subsystems  of 
the  ■ nstrumeni , a test  was  ecmceived  to  demonstrate  the  overall  capabilities 
of  the  frequency  tracking  function.  An  artificial  whistler  generator  was 
constructed  by  connecting  a discharging  R-C  circuit  to  a voltage- controlled 
oscillator  to  produce  a sine-wave  sweeping  from  10  kHa  down  to  zero 
frequency  with  an  exponential  profile  resembling  a whistler.  The  output  of 
this  oscillator  was  mixed  with  a white-noise  generator  to  simulate  VLF  back- 
ground noise.  By  restricting  the  noise  bandwidth  to  10  kHz  with  a low -pass 
filter,  a quantitative  measure  of  the  noise  power  spectral  density  could  be 


made. 


With  the  whistler-and-noise  simulator  connected  to  the  instrument  and 


the  controls  adjusted  so  that  the  receiver  would  intercept  the  whistler  at 
9 . 2 kHz  and  track  it  down  to  500  Hz , the  speed  of  the  whistler  and  the  level 
of  the  noise  were  increased  until  the  limits  of  the  tracking  capability  were 
reached . The  initial  frequency  excursion  rate  of  the  whistler  was  32  kHz  per 
second  and  the  noise  power  spectral  density  was  -32  dBm/Hz  for  this  condition. 
(The  signal  level  was  0 dBm . ) This  corresponds  to  a signal -to-noise  ratio 
of  “8  dB  in  the  full  10  kHz  VLF  bandwidth , or  +6 . 7 dB  in  the  340  Hz  IF 
bandv/idth . Fig.  4.1  shows  analog  chart  records  produced  during  this 
experiment  demonstrating  the  tracking  action  and  the  <x»rrect  measurement  of 
the  signal-to-noise  ratio  by  the  parameter  extraction  circuitry . 

Further  experimentation  revealed  that  higher  tracking  rates  and  lower 
signal-to-noise  ratios  could  be  tolerated  once  the  instrument  was  locked  onto 
a signal . approaclilng  the  theoretical  limits  derived  in  Section  II . For 
practical  purposes , however . when  an  unknown  signal  must  be  acquired  by 
the  instrument,  the  limits  determined  by  the  test  are  more  realistic. 

To  teat  the  ability  of  the  instrument  to  track  whistlers  and  other  VLF 
emissions  in  an  actual  spectral  environment . tape  recordings  from  Vandenberg, 
CalifoKiia;  Eights  Stalltm,  Antarctica;  and  Rbberval,  Quebec  wore  played  as 
inpvrt  to  the  equipment . Successful  tracking  was  achieved  <m  strong  whistlers , 
weak  whistlers,  nose  whistlers , and  roagnetoapherlcally  propagated  signals 
from  the  transmitter  at  Slple  Station.  Antarctica  received  at  Eoberval.  Since 
a demonstration  of  the  frequency-tracking  capability  is  even  more 
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LABORATORY  TRACKING  TEST 
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FIGURE  4.1  LABORATORY  TRACKING  TEST . (a)  Frequency  of 

electronically  simulated  whistler.  Initial  slt^e  is  -32  kHz 
per  second,  (b)  Frequency  of  tracking  receiver,  showing  nature  of  the 
'“whistler”  « 9.2  kHz  and  successful  tracking  down  to  lower  operating 
frequency  limit,  (c)  Signal  to  noise  ratio  indicated  by  parameter  extraction 
circuitry . The  indication  agrees  well  with  the  measured  signal  to  noise 
ratio  of  6 . 7 dB . 
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convincingly  presented  ^om  actual  real  time  field  operations , which  are 
discussed  in  the  following  sections , the  data  records  from  these  laboratory 
experiments  will  not  be  reproduced  here . 

2.  Direction-Finding . 

The  Instrumental  accuracy  of  the  direction-finder  was  measured 
using  simulated  directional  signal  inputs . Inverted  and  noninverted  versions 
of  a sine-wave  signal  were  generated  by  operatiomd  amplifiers . The  vertical 
anteima  (channel  1)  input  was  connected  to  the  noninverted  signal . By 
connecting  the  two  loop  inputs  (channels  2 and  3)  either  to  the  noninverted 
signal , the  inverted  signal , or  ground . signal  directions  of  arrival  of  0 , 
t45°,  t90°»  ±135*^  and  180°  could  be  simulsied.  The  instrumental  errors  for 
these  angles  > measured  over  a 40  dB  range  of  signal  levels , are  shown  in 
Table  4.1.  Overall  the  RMS  error  for  all  the  measurements  was  1 . 2S° . 


TABLE  4. 1 mSTRUMENTAL  ERRORS  OF  DIRECTION-FINDER. 


Simulated  Angle 

-136° 

-00° 

-46° 

0° 

+45 

+00° 

+135° 

180“1 

Relative  Signal 
Level 
(dB) 

Instrumealal  Error  (D^reos) 

_ _ j 

0 

0.0 

+0.6 

+0.5 

-0. 6. 

+1.0 

0.0 

0.0 

+1.0 

-10 

+0.5 

'-0.5 

-1.6 

-0.6 

+S.0 

-0.5 

-1.5 

0.0 

-20 

-0.6 

-1.0 

-0.6 

-l.O 

+1.0 

-1.0 

-0.6 

-0.6 

-30 

-2.8 

-l.S 

o.-o 

O 

1 

-1.5 

-l.u 

+0.5 

+1.0 

-40 

+0.8 

+4.0 

+0.5 

-1.5 

+1.0 

-2.0 

-2.0 

-2.0 

Again , as  in  the  case  of  the  freqlfency-trr  iking  capability , the  most 

convincing  demonstration  of  the  direction-finding  operation  of  the  instrument 

«** 

is  provided  by  the  actual  real-time  records  presented  in  the  sections  which 
follow . 

B.  FIELD  OBSERVATIONS  AT  STANFORD,  CALIFORNIA 

The  prototype  frequency-tracking  direction-finder  was  installed  at  the 
Stanford  VLF  observiiig  site  and  operated  during  the  period  from  28  March  to 
9 June  1975.  Observations  were  made  on  subionospherically  propagated  VLF 
transmitter  signals  and  on  whistlers . 

1.  Site  Corrections 

As  no  survey  data  were  available  for  the  antenna  system  at  Staitford , 
no  corrections  could  be  applied  for  antenna  misalignment . Earlier  work  at 
this  site  by  Cousins  [l972]  indicated  site  errors  up  to  5°  according  to 
bearings  taken  on  VLF  stations.  The  position  of  the  feedlinos  running  from 
tlie  antennas  to  the  equipment  building  had  since  been  modified,  however,  so 
no  attempt  was  made  to  apply  his  data  as  a correction  factor. 

2-  VLF  Transmitters 

On  19  March  1975.  bearings  were  taken  on  four  VLF  transmitters. 

Portions  of  the  chart  record  showing  the  bearing  indications  for  the  stations 
Omega/Hawaii,  Omega/Ncrth  Dakota,  NPG,  and  NAA  are  presented  in  Fig.  4.2. 
The  bearing  angles  extracted  from  the  chart  are  wmpared  with  the  calculated 
great  c-ircle  beatdngs  of  these  stmions  in  Table  4.2.  The  errors  may  be 
attributed  mainly  to  situ  errors  and  antentta  misalignment . and  their  magnitudes 
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FIGUES  4.2  DIRECTION-FINDING  ON  VLF  STATIONS  RECEIVED  AT 
STANFORD . Chart  data  show  direction  of  arrival  and 
relative  signal  strength . (a)  OMEGA  atfations  in  Hawaii  and  North  Dakota 

Iransraitting  in  sequence  on  lx. 33  kHz.  (b)  NPG,  18.6  kHz.  (c)  NAA, 
17 . 3 kHz . Results  are  given  in  Table  4.2. 


TABLE  4. 2 DIBECTION'FlHDiNG  RESULTS  ON  \nLF  STATIONS  HSCEiVED 
AT  STANFORD, 


Station 

Frequency 

L-ocatlon 

indicated 

Bearing 

True 

Bearing 

1 

Error 

OMEGA/H. 

li.33 

Hawaii 

257.0® 

262. 6® 

+4.4® 

OMEGA/N.D. 

11.33 

N.  Dakota 

59.4® 

55.6® 

+8.8® 

NPG 

IS.  60 

Washington 

1.7® 

1.0° 

+0.7® 

NAA 

17.80 

Maine 

88.4® 

62.0® 

+6.4® 

are  close  to  those  of  Cousins'  results . For  the  signals  from  NAA  and  Omega/ 
North  Dakota,  polarization  error  may  also  have  been  a factor,  as  sunset  was 
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to 


occurring  along  the  path  and  this  condition  was  shown  by  Horner 
be  conducive  to  such  errors . 

To  further  demonstrate  the  usefulness  of  the  instrument  in  direction- 
finding  studies  of  polarization  error  on  the  signals  from  VLF  transmitters 
continuous  records  of  the  bearings  for  several  stations  were  made  over 
24-hour  periods . In  agreement  with  Homer's  results , the  daytime  errors 
were  low , rising  at  sunset  and  varying  cyclically  through  the  night , and 
falling  off  at  sunrise,  A particularly  interesting  example  for  station  NPG 
showing  a drainatic  change  in  the  behavior  of  the  polarization  error  just 
before  sunrise  is  presented  in  Fig.  4.3. 


FIGURE  4.3  POLARIZATION  ERROR  BEHAVIOR  IN  THE  HOURS  BEFORE  LOCAL 
SUNRISE  (1246  UT)  FOR  STATION  NPG  RECEIVED  AT  STANFORD. 


3.  Whistlers 


The  VLF  observing  site  at  Stanford  was  unattended  during  the 
whistler  monitoring  program.  Human  attention  was  limited  to  one  brief  visit 
every  two  days  < at  which  time  the  operator  could  mount  one  new  reel  of  tape 
(holding  U hours  of  data)  and  change  the  equipment  adjustthents  or  select 
a diffex^nt  time  period  for  the  next  monitoring  run . The  actual  recordhtg  was 
triggered  automatiealiy  by  a clock  to  occur  during  an  early  morning  hour, 
sfnce  that  period  was  known  to  be  favorable  for  whistler  reception . These 
conditions  made  it  imperative  that  the  operator  have  immediate  data  on  the 
rpsults  of  the  previous  run  so  that  he  could  intelligently  modify  the  equipment 
adjustments  or  the  recording  time  period  for  the  next  run. 

To  fulfill  tiiis  requirement,  the  signal  presence  detection  capability  of 
the  instrument  was  applied  to  cause  a pen  deflection  on  a chart  recorder 
ramting  at  a very  slow  speed.  With  this  arrangement , the  operator  could 
scan  an  entire  two-day  period  quickly  to  determine  if  signals  were  being 
traekad . if  any  occurred  during  the  li  hour  tape  recording  period , they 
could  be  located  immediately  and  aurally  reviewed  without  listening  to  or 
spectrum  analy  sing  the  entire  U pe. 

A total  of  18  recordings  were  made , containing  27  hours  of  VLF  data. 
Eleven  whistlers  were  found  which  ’vex-e  tracked  by  the  equipment , most  of 
them  in  the  latter  part  of  the  progrwi  when  the  e';Uipment  settings  had  been 
optimized.  Two  examples  are  shown  in  Fig.  4.4.  While  both  bearing 
indications  show  cousiderahle  variebiiity , the  example  of  9 May  shows  a 


100 


9 MAY  75  II  MAY  75 
llll;37UT  1151:54  UT 


FIGURE  4.4  WHISTLERS  RECEIVED  AT  STANFORD,  (a)  VLF  spectrum 
showing  whistler  on  9 May  1975 . (b)  Spectrum  analysis  of 

VLF  mixed  with  FM  data  channel  showing  tracking  of  whistler  frequency . 

(c)  Signal  strength  and  direction  of  arrival  data  for  the  whistler . Bearing 
data  show  considerable  variability  but  a predominance  toward  northeast  to 
east,  (d) , (e) , and  (f)  Similar  records  for  a whistler  on  11  May  75,  showing 
predominantly  western  bearing.  Both  examples  show  evidence  of  multiple 
components  which  may  be  responsible  for  varfability  of  bearing  indication . 
(Blanking  of  DF  trace  in  the  absence  of  a signal  was  not  included  at  this  time . ) 


general  easterly  direction » while  the  whistler  of  11  May  displ^s  generally 
west  to  southWeat  bearings . The  vatiability  within  each  example  is 
attributable  to  multiple  Whisfler  components , which  are  clearly  demonstrated 
by  the  rather  diffuse  spectral  trace  of  the  whistlers . The  performance  of 
the  direction-finder  in  the  presence  of  multiple’  signals  was  discussed  in 
Section  II.C.2,  where  it  was  shown  that  the  average  bearing  would  tend  to 
be  that  of  the  strongest  signal  component , with  periodic  variations  about  the 
average  having  a magnitude  dependent  on  the  relative  amplitude  and  bearing 
angles  of  the  interfering  components , and  with  a periodic  rate  corresponding 
to  the  difference  frequency  between  the  primary  component  and  the  interference 
While  speciffc  components  cannot  be  identiffed  in  the  examples , the  bearing 
displays  are  consistent  with  such  an  explanation . 

C . FIELD  OBSERVATIONS  AT  ROBERVAL , QUEBEC 

The  prototype  frequency-tracking  direction-finder  was  installed  at 
Roberval , Quebec  and  operated  during  the  period  from  June  until  the  end 
of  August  1975.  Observations  were  made  on  VLF  transmitters , natural  VLF 
activity  including  whistlers , emissions , and  chorus , and  on  artificially 
stimulated  emissions  produced  by  the  signals  from  the  Stanford-operated 
transmitter  at  Siple  Station,  Antarctica. 
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X.  Site  Corrections 


A survey  of  the  mounting  points  for  the  loop  antennas  revealed 
their  orientation  to  be  toward  magnetic  north  rather  than  true  north  . 
Accordingly*  the  analysis  of  misalignment  effects  presented  in  Secdion  II. C. 4 
was  applied  based  on  the  survey  data.  The  quantities  6 = -19.92^  and 
* = -17.56°  were  obtained  from  the  survey.  From  equation  (2.88) , the 

A 

measured  direotion-of-arrival  6 . contains  an  error. 


A8  = 18.74°  + 1.18°  cos  20. 


* 

Solving  for  6 In  terms  of  6 . 


e = 0 - 18.74°  - l.l8°coa2  (9-  18.74°). 

The  raw  data  records  presented  herein  show  tlie  indicated  bearing  in  all 
tabulated  data  and  discussions  of  results,  however,  the  corrected  bearings 
are  used. 

2*  VLF  Transmitters 

Bearing  data  were  taken  on  stations  NPQ , QBR,  and  four  Omega 
stations  which  could  be  received  at  Roberval.  The  chaxt  records  produced 
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during  these  observations  are  shown  in  Fig.  4.5 » with  the  corrected  bearings 
and  a comparison  to  the  calcul^ed  great  circle  bearings  given  in  Table  4,3. 
Bearings  on  the  VLF  stations  were  found  to  be  quite  accurate  and  highly 
repeatable  from  day  to  day . 


TABLE  4. 3 DIRECTION-FINDING  RESULTS  ON  VLF  NATIONS  RECEIVED 
AT  ROBERVAL,  QUEBEC. 


Station 

Frequency 

Location 

Indicated 

Bering 

Corrected 

Bearing 

True 

Bearing 

Error 

GBR 

16.0 

England 

74.3® 

68.0® 

56.8® 

-0.8® 

NPG 

18.6 

Washington 

307. S® 

289. 7® 

288.7® 

+1.0® 

OMEGA/NOR 

10.2 

Norway 

58.4® 

34.3° 

34.3® 

0.0® 

OMEGA/TRI 

10.2 

Trinidad 

189.7® 

169.0® 

163.0® 

+6.8® 

OMEGA/H. 

10.2 

Hawaii 

297,9® 

280.3® 

281.1® 

-0.8® 

OMEQA/N.D. 

10.2 

N.  Dakota 

286.2® 

268.6® 

273.1® 

-4.5® 

3,  Whiatlers  and  Other  N 
Conditions  for  observations  at  Roberval  were  considerably  more 
favorable  than  at  Stanford.  A greater  abundance  of  VLF  activity  was 
available » and  the  site  was  constantly  manned  during  regular  daily  operating 
periods , permitting  the  instrument  to  be  adjusted  to  meet  the  chmiging  signal 
environment.  In  additioii,  a general  confirmation  of  direction'finding  results 
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FIGURE  4 . 5 DIRECTION-FINDING  ON  VLF  STATIONS  RECEIVED  AT 

ROBERVAL . Chart  data  show  direction  of  arrival , measured 
signal -to-noise  ratio « and  relative  signal  strength,  (a)  Station  GBR  on 
16 .0  kHz . Cb)  NPG  ,18.6  kHz . (c)  OMEGA  stations  in  Norway . Trinidad , 

Hawaii,  and  North  Dakota  transmitting  in  sequence  on  10.2  kHz.  Results  are 
given  in  Table  4.3. 


was  made  possible  by  comparing  the  bearings  with  the  indicated  magnetic 
latitude  ''f  the  signal  paths  based  on  whistler  nose  frequency  analysis . For 
the  purposes  of  this  comparison , a propagation  model  set  up  by  Seely  [iSTs] 
predicting  nose  frequencies  of  3.5  to  3.6  KHz  at  the  Roberval  latitude  was 
employed.  Whistlers  with  higher  nose  frequencies  would  be  expected  to 
show  southerly  bearings , while  those  with  lower  nose  frequencies  should 
be  found  to  the  north.  Signals  with  comparable  nose  frequencies  mig^t 
display  any  bearing,  although  if  the  bearing  were  near  north  or  south  the 
signal  would  have  to  be  nearly  overhead:  in  this  case  other  signs  mi^t  be 
present . such  as  high  polarization  error,  indicating  a high  elevation  angle 
of  arrival . For  non-whistler  signals  believed  to  propagate  in  the  same  ducts 
as  whistlers , such  as  triggered  emissions,  the  nose  frequencies  can  be  taken 
from  whistlers  occurring  in  the  same  time  period  for  the  purposes  of  this 
comparison . 

An  example  of  nose  whistlers  and  associsted  emissions  is  given  in 
Fig.  4.6.  In  many  cases . the  instrument  was  surprisingly  successful  in 
tracking  the  frequency  of  weak  whistler  components  and  emissions;  this  is 
demonstrated  by  the  mixed  record  (b)  in  which  components  which  were  faint 
in  the  VLF  spectnun  are  made  clesily  visible  by  the  superposition  of  the 
tracker  frequent^  trace . Within  this  record  one  can  see  oases  of  highly 
variable  bearings  (1.0  second)  and  relative^  consistent  bearings  (S.5 
seconds) . These  differences  are  attributable  to  multiple-  veraxis  single- 
component reception . Later  examples  will  illiuninate  this  more  clearly. 
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FIGURE  4.6  NOSE  WHISTLERS  AND  ASSOCIATED  EMISSIONS . U>  VLF 

spectrum  analysis  showing  whistlers  and  associated  emissions 
(b)  Spectrum  analysis  of  VLF  mixed  with  tracking  frequency  trace 
demonstrating  tracking  o?  signals . <o)  Sl«mal  strength  and  direction-of- 

arrival  data.  The  signals  display  a predominantly  south-southwest  bearing, 
in  agreement  with  the  indication  by  the  whistler  nose  frequencies  (''•4.3  kBe) 
that  the  paths  are  southward  of  Roberval 


A 


The  consistent  bearings  indicate  a south-southwest  (corrected) direction, 
in  agreement  with  the  nose  frequencies  at  4.3  kHz  which  indicate  an  exit 
point  south  of  Reberval . 

In  the  next  example.  Fig.  4.7,  taken  one  minute  later  on  the  same  day , 
a different  perspective  is  provided.  A 30-second  time  period  is  displayed  to 
demonstrate  the  consistent  south-southwest  (corrected)  bearings  indicated 
for  repeated  excitations  of  the  same  whistler  ducts . Also  visible  on  this 
record  are  some  strong  spheric  pairs  (at  7 seconds  and  13  seconds)  which 
caused  the  two-hop  whistler  pairs  at  11  seconds  and  17  se<»nds . Bearings 
are  indicated  for  these  spherics  also:  south  and  then  northwest  &r  the 
first  pair;  south  and  then  northeast  for  the  second  pair.  The  very  strong 
whistlexvassociated  emission  at  27  seconds  is  the  subject  of  the  next  example. 

The  VLF  spectrum  in  Fig.  4.8(a)  shows  some  weak  whistler  comiKments 
and  an  intense  emission  triggered  by  one  of  these . The  intensity  of  the 
emission  is  also  evident  the  signal  strength  trace  in  (c) . In  accordance 
with  the  high  signal-to-noise  ratio  which  exists  for  such  a strong  signal , the 
bearing  indication  is  extremely  steady . The  direction  Is  identical  to  the 
whistlers  of  Pig.  4.7,  evidence  that  the  emission  has  propagated  through 
the  same  duct  as  the  whistlers . 

Several  cates  were  observed  in  which  the  omissions  associated  with 
whistlers  displayed  better  coherency  than  the  whistlers  themselves . An 
illustration  of  this  obser^^atlon  is  presented  in  Fig.  4.9.  The  spectrum 
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FiOURB  4 . 7 MULTIPLE  WHISTLERS  AHD  ASSOCIATED  EMISSIONS . 

(a)  VLF  spectrum  analysis  showing  repeated  whistlers . The 
causative  spherics  are  also  visible . (b)  Mixed  VLF  and  frequency  tracker 

idtowing  tracking  of  whistle» . (c)  Signal  strength  and  <Urection-of- arrival 

data.  Bearings  are  consistent  between  the  repeated  whistler  excitations. 
The  bearing  of  the  causative  spherics  are  also  displayed. 
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FIOUHE  4.8  WmSTLERS  AND  ASSOCIATED  EMISSION . (a)  VLF  spe<Jtrum 

analysis . (b)  Mixed  VLF  spectrum  aitd  frequency  tracking 

trace  showing  instrurnsnt  tracking  emission  and  weaker  whistler  components 
(c)  Signal  strength  trace  demonstrates  high  intensity . Bearing  indication  is 
extremely  steady  in  accord  with  high  signal-to-noise  ratio  existing. 
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FIQUiiE  4.9  A DIFFUSE  WHISTLER  WITH  STRONG  PRECURSOR  EMISSION . 

(»)  Spectrum  analysis  showing  diffuseneas  of  whistler. 

(b)  Mixed  VL?  and  tracking  frequency  record.  Diffuse  nature  of  whistler 
causes  ♦racker  to  follow  oniy  the  general  energy  maximum . (c)  Bearing 

indications  showing  gradually  changing  bearings  from  southwest  to  south- 
east on  precursor  and  extremely  confused  bearing  on  whistler. 


analysis  in  (a)  shows  a strong  precursor  emission  and  an  extremely  diffuse 
multicomponent  whistler.  The  whistler  was  so  diffuse  that  the  frequency 
tracker  could  only  follow  the  general  energy  maximum  of  the  spectrum;  the 
bearing  indication  in  (c)  is  accordingly  very  confused.  The  precursor 
displays  better  coherency  . Its  bearing  indication  shows  less  rapid  variation, 
although  a pronounced  variation  with  frequency  is  evident  as  the  bearing 
drifts  from  southwest  > through  south , to  a final  reading  near  southeast . This 
frequency  variation  may  be  considered  evidence  of  polarization  error,  which 
in  turn  suggest;;  a moderately  high  elevation  angle  of  arrival . As  a broader 
data  base  is  gathered,  direction-finding  studies  may  prove  of  value  in 
determining  the  origin  end  propagation  path  of  the  less  understood  emissions 
such  as  precursors . 

VLF  emissions  not  associated  with  whistlers  were  also  observed  at 
Roberval.  On  10  July  1975 , a buildup  of  chorus  activity  occurred.  Initially, 
the  instrument  had  been  set  up  to  monitor  for  whijtlers  at  3.3  kHz,  In 
Fig.  4.1*',  two  hook-type  emissions  associated  with  chorus  are  tracked,  as 
shown  by  the  tracking  frequency  trace  in  (b) , at  0 seconds  and  4.5  seconds. 
Bearing  indications  in  (c)  are  steady  and  consistent , displaying  a (corrected) 
south -southwest  direction.  The  operator  then  retuned  the  instrument  to 
monitor  at  2 kHz , causing  tracking  to  occur  on  the  chorus  elements  themselves 
as  shown  in  Pig.  4.11.  The  bearing  indications  changed  to  west -southwest 
and  displayed  somewhat  more  variability , as  one  might  expect  from  the 
multiplicity  of  components  present  in  the  chorus  activity . Taken  approximately 
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FIGURE  4 . 10  HOOK-TYPE  EMISSIONS  ASSOCIATED  WITH  CHORUS  . 

(a)  Spectrum  analysis  showing  chorus  and  associated  hooks 
at  0 and  4.5  seconds . (b)  Mixed  VLF  and  tracker  record  demonstrating 

capture  of  hooks  by  frequency -tracker  waiting  at  3 . 3 kHz . (d)  Bearing 

indications  showing  steady  south-southwest  (corrected)  direction  for  the 
hooks. 
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twelve  minutes  later,  the  example  shown  in  Fig.  4.12  shows  a hemming  trace 
with  a considerable  amount  of  both  re^id  and  slow  variations  indicative  of 
polarization  error.  A possible  explanation  is  a drift  of  the  ducts  carrying  the 
chorus  to  a nearly  overhead  location . Drifts  of  this  nature  are  of  interest 
in  stupes  of  magnatospheric  convection  inentioned  in  Section  I.B.I  as  a 
motivating  factor  for  direction-finding  . 

A more  detailed  example  of  temporal  vmdations  in  the  apparent  angle 
of  arrival  is  provided  by  Fig.  4.13,  in  which  three  successive  40-second 
samples  t^en  3 minutes  apart  on  28  July  are  displayed.  In  order  to  display 
the  signals  and  directional  data  compactly , a single  0 to  20  kHz  spectrum 
analysis  of  the  VIF  data  mixed  with  the  FM  data  channel  is  utilized.  Sample 
(a)  at  1539  hours  shows  whistlers  and  emissions  displaying  southerly  bearings . 
In  sample  (b) , the  bearing  data  are  frequently  broken  up  due  to  the 
instrumental  scale  discontinuity  at  due  south.  At  1545  hours , in  sample  (c> , 
the  signals  have  changed  character,  displaying  more  diffuseness  and  a lower 
whistler  nose  frequency;  simultaneously , the  bearing  indications  have  moved 
to  a south-southeast  direction.  Recalling  horn  Figs  .4.6  and  4.7  the  south- 

^A  scale  discontinuity  is  unavcddable  in  representing  a continuous  angular 
variable  on  a finite  linear  scale . For  signals  incident  near  the  discontinuity 
the  indication  may  jump  erratically  between  the  extremes  of  the  linear  scale. 

If  the  choice  of  south  as  the  discontinuity  proves  inconvenient , the  phase  of 
the  vertical  antenna  signal  can  be  reveraedto  reflect  all  angles  through  the 
origin,  placing  the  discontinuity  at  north  instead. 
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FIGURE  4.12  CHORUS,  (a)  VLF  spectrum  analysis,  (b)  Mixed  record 

showing  tracking  of  chorus  elements . (c)  Bearing  indications 

show  both  slow  and  rapid  variations  suggestive  of  polarisation  error  and  a 
hi{^  elevation  angle  of  arrival . Compared  with  Fig.  4.11,  a drift  from  west 
to  overhead  may  be  Indicated, 
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southwesit  besring's  at  1530  hours , a ;^ssibla  picture  of  a nov&koast  drift  is 
evident . While  farther  analysis  and  a broader  data  base  are  required  before 
direct  conUrmatiGn  of  duct  drift  th^kes  can  be  pres^tedr^^b  examples  do 
provide  evidence  of  the  usefulness  of  the  direction-Snder^  -sbeh^^to^ 
tions. 

4*  Artificially  Stimulated  ghdasions  fet^^Stele  jgtation , Antarctica. 

Roberval  is  located  at  the  magnetically  coi^di^ate  point  to  Siple 
Station,  Antarctica,  where  Stanford  University  operates  a 100  kW  VLF 
transmitter  for  i^robing  of  the  magnetosphere  [Helbwell  and  Katsufrakis , 

1974j  . Signals  transmitted  at  Siple  stimulate  emissions  and  display  amplifica- 
tion when  received  at  Roberval . The  operation  of  the  prototype  frequency- 
tracking  direction-finder  at  Roberval  coincided  with  a campalj^of  experimenta- 
tion employing  the  Siple  transmitter , and  signals  from  Siple  were  observed  on 
several  occasions . The  programs  transmitted  included  CW  pulses , frequency 
ramps , and  frequency  sawtooth  sequences . 

A typical  reception  of  frequency  ramps  from  Siple  is  illustrated  in 
Fig.  4.14.  In  (a) , the  VLF  spectinun  analysis  shows  the  repeated  ramps 
followed  by  a descending  "staiTCase”  sequence.  Also  visible  are  the  idler 
pulses  used  during  key-up  conditions  to  maintain  the  load  on  the  transmitter 
but  which  are  amplified  much  less  than  the  ramps.  The  mixed  VLF  and 
tracker  record  (b)  demonstrates  the  successful  tracking  of  the  signal 
frequency  by  the  instrument , including  the  ramps  and  descending  staircase . 
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The  bearing  trace  in  (c)  shows  a characteristic  frequency-dependent 
variation  during  the  first  part  of  each  pulse  which  may  indicate  polarization 
error  due  to  a relatively  high  angle  of  arrival . The  bearing  signature  repeats 
fbr  each  successive  pulse.  The  bearings  on  the  staircase  steps  agree  with 
the  bearing  for  the  corresponding  frequency  within  the  ramps , although  the 


changes  are  harder  to  detect  due  to  the  noiaineae  of  the  trace . From  the 


north  to  northwest  bearing  and  the  presence  of  polarization  error,  an  exit 
point  near  the  station  and  to  the  north  would  be  inferred.  This  is  consistent 
with  analyses  of  nose  whistlers  from  the  same  period,  which  show  nose 
frequencies  of  approximately  3.2  kHz. 

Additional  examples  from  the  same  date  are  provided  in  Figs.  4.15 
and  4.16.  In  4 . 15 , an  apparent  short  term  change  in  the  polarization  error 
characteristics  occurs . At  1108: 51,  the  whistler  and  Siplo  signal  show  the 
same  frequency-dependent  bearing  variations  as  in  Fig.  4.14.  At  1109: 36, 
however,  the  bearing  variation  with  frequency  has  changed  in  character. 

This  rapid  change  is  consistent  with  the  sensitivity  of  the  polarization  ellipse 
to  small  changes  in  elevation  angle  as  the  multiple  rays  from  the  exit  point  to 
the  receiver  shift  relative  to  one  another.  The  example  of  Fig.  4.16 
illustrates  the  effect  of  multipath  versus  single  component  reception  conditions 
on  the  bearing  indication . In  the  example  at  1109: 27 , the  spectrum  analysis 
at  (a)  displays  multipath  conditions  on  the  Siple  ramp  and  nose  whistler. 

The  tracking  frequency  trace  in  (b)  shows  a characteristic  fuzziness  in- 
dicating low  coherence . At  (c)  the  directional  indication , although  indicating 
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FIGURE  4 . 15  SIPLE  SIGNALS  SHOWING  CHANGING  POLARIZATION  ERROR 
(a)  and  (c)  Spectrum  analysis  of  VLF  mixed  wiGi  tracking 
frequency  indication  for  two  ramps  separated  by  approximately  <me  minute 
(b)  and  (d)  Direction-of-arrival  data  showing  a distinct  difference  in  the 
bearing  behavior  through  the  ramp  duration.  A change  in  polariaidion 
conditions  may  be  inferred. 


L2 


R0  7AU6  75 

ll09c27UT  flf4;22UT 


FIGURE  4. 16  SIPLE  SIGNALS  SHOWING  EFFECTS  OF  MULTIPATH . 

(a)  Specarum  analysis  of  Siple  ramp  and  a nose  whistler 
displaying  multipath  conditions,  (b)  Tracking- frequency  trace  and  VLF 
mixed  record  showing  successful  tracking  but  noisy  trace  indicative  of  low 
coherence . (c)  Bearing  indications  show  a consistent  average  but 

considerable  rapid  variation.  The  example  of  1114  (d)  shows  less  multipath 

(e)  Frequency  tracking  trace  for  latter  part  of  ramp  is  sharply  defined. 

(f)  Bearing  indication  is  more  consistent . 
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a consistent  average > is  subject  to  a great  deal  of  rapid  periodic  variation. 

In  contrast,  the  example  at  1114:22  shows  less  multipath  in  the  spectrum 
analysis  (d) , and  a more  clearly  defined  tracking  frequency  trace  in  (e) 
for  the  lattem  portion  of  the  Siple  ramp . The  corresponding  bearing 
indio^ion  (f)  shows  less  variability . 

An  even  more  extreme  case  of  multipath  is  given  in  Fig.  4.17.  The 
bearing  indlce^on  (c)  is  extremely  confused  except  for  brief  instants  at 
2.0  seconds  and  12.0  seconds.  A study  of  the  spectrum  (a)  and  the  tracking 
frequency  record  (b)  shows  that  relatively  coherent  triggered  emissions 
carried  the  tracker  foequency  away  from  the  multicomponent  sawtooth  signal 
from  Siple  at  these  times . The  tracking  frequency  trace  as  recorded  on  a 
chert  recorder  (b)  shows  a striking  change  from  a fuzzy  character  to  a clean 
trace  during  these  emissions . indicating  the  higher  degree  of  frequency 
coherent^  for  the  emissions . 

A giimpae  at  another  application  the  instrument  is  provided  in 

Fig  .4.18.  In  addition  to  the  familiar  quantities  of  direction-c^-arrival  (c) , 

signal-to-nolse  ratio  (d),  signal  strength  (e),  and  tracking  frequency  (f). 

the  real  and  imaginary  parts  of  the  r«4to  (-H  /H  ) (which  are  also  measured 

X y 

by  the  instnunent)  are  shown  at  <b)  and  (a).  Noting  that  the  imaginary 

part  ia  negative  and  greater  than  the  real  pert  during  reception  of  the  ramps 

frcmi  Siple , one  can  conclude  that  the  polarization  is  essentially  circular.  The 

handedness  can  also  be  ascertained  from  the  sign  of  U /U  : as  ('H  /H  ) is 

X y X y 


FiaURE  4 . 17  SAWTOOTH  SIGNAL  FROM  SIPL2  AND  TRIGGERED  EMISSIONS . 

<s)  VLF  spectrum  analysis  showing  multipath  reception  of 
Slplo  sawtooth,  (b)  Tracking  frequency  chart  record  displays  triggered 
emissions  trwiked  at  2.0  and  12.0  seconds.  Trace  character  becomes  sharper 
indicating  higher  coherence  ter  these  emissions . (c)  Bearing  record  is 
extremely  confused  due  to  multipath , with  the  exception  of  the  two  emissions . 
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FIGURE  4.18  CH/.RT  RECORD  INDICATING  POLARIZATION  OF  RECEIVED 
SIPLE  RA.MPS.  (a)  and  (b)  Measurements  of  the  real  and 
imaginary  parts  of  the  ratio  (-H  /H  ) are  recorded,  (c)  Bearing 
signature  shows  C'^nsiderable  frequ'liicy  variation  indicative  of  high  elevation 
angle,  (f;  Frequency  record  showing  tracking  of  i-amps.  The  polarization 
of  these  signals  is  essentially  right-hmid  circular  as  deduced  from  (a)  and  (b) 


indicated  as  negative,  H /H  must  be  positive;  therefore,  H leads  H . 

The  polarization  must  be  right  hand  circular  Gooldng  down  at  the  receiving 
site , in  the  direction  of  propagation) . This  agrees  with  Helliwell's  [l965] 
description  of  the  polarization  of  whistler-mode  waves . The  indication  of 
circular  polarization  strongly  suggests  a high  elevation  angle  of  arrival  and 
a near-overhead  duct . Again , this  is  in  agreement  with  hose  frequency 
observations  of  3.1  to  3.6  kHz  on  this  date. 
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V.  CONCLUSIONS 


A.  SUMMARY  OF  RESULTS 

A theory  and  technique  fo?  tracking  the  frequency  of  whistlers  and 
related  VLF  emissions  and  extracting  direction-of-arrivaJ,  data  in  real  time 
has  been  presented.  A prototype  instrument  embodying  these  innovations 
has  been  designed  and  constructed , and  tested  both  in  the  laboratory  and 
at  field  observing  sites . 

Laboreiory  tests  of  the  tracking  capability  show  that , with  a 340  Hz 
IF  filter  bandwidth , signals  changing  frequency  at  up  to  32  kHz/second  can 
be  captured  and  tracked  at  signal-to-noise  ratios  as  low  as  6 . 7 dB . At 
field  observing  sites , successful  tracking  was  achieved  for  whistlers , both 
above  and  below  the  nose  frequency , and  on  a variety  of  emissions , chorus , 
and  Siple  transmitter  signal  formats , in  the  presence  of  both  manmade  and 
atmospheric  noise . 

Laboratory  tests  of  the  direction- finding  capability  show  an 
instrumental  accuracy  of  1 . 25^  RMS . Field  observations  on  several  VLF 
transmitting  stations  show  repeatable , accurate  results , with  the  average 
magnitude  of  the  error  being  2.9°,  most  of  which  is  attributable  to  site 
error. 

Direction-finding  bearings  were  obtained  on  whistlers , sphex’ics , 
chorus , Siple  signals , and  related  emissions . The  behavior  of  the  bearing 
indications  agrees  with  the  theoretical  prediction  of  system  behavior  in  the 
presence  of  noise , multipath , and  polarization  error.  Under  favorable 
conditions,  repeated,  consistent  bearings  were  obtained. 


The  bearing  results  on  those  signals  were  omsistent  with  the  exit 
point  latitudes  predicted  from  whistler  nose  frequency  analysis , without 
the  necessity  of  invoking  deviations  of  the  ray  path  from  the  field  line  of  the 
duct. 

Some  of  the  data  display  a temporal  change  in  the  indicated  bearing 
which  may  be  indicative  of  a duct  drift. 


B.  CONCLUSIONS 

The  new  technique  for  VLF  direction-finding  presented  herein  fills 
a gap  in  the  range  of  techniques  currently  available  or  under  development . 
The  directional  data  are  of  higher  accuracy  and  temporal  resolution  than 
those  produced  by  the  goniometer  technique . The  instrument  is  of  use  against 
a wider  variety  of  signals  and  wave  incidence  conditions  than  the  multi- 
channel techniques  of  Cousins  [l972  or  Tsuruda  and  Hayashi  [l974]  . 

Finally , the  availability  of  the  information  in  real  time  and  with  a minimum 
of  human  analytical  effort  represents  an  advajiitage  over  both  of  the  other 
methods . The  practicality  of  the  instrument  has  been  demonstrated  by  the 
success  of  the  field  observation  program . The  frequency-tracking  direction- 
finder promises  to  yield  data  hitherto  unavailable  on  the  location  of  the  ducts 
which  carry  whistlers  and  related  signals. 
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C.  FUTUKE  APPUCATIONS 

A number  of  different  opportunities  application  of  the  instrument 
are  apparent . To  further  exploit  the  direction-finding  capability , it  is 
recommended  that  a multistation  network  be  employed . Two  stations  would 
permit  location  of  sources  by  triangulation . Three  of  more  stations  would 
provide  a further  check  on  the  system  accuracy  by  examining  the  eon* 
vergence  of  the  bearings  to  a single  point . 

The  capability  of  the  instrument  for  measuring  the  apparent  polarizatlDn 
of  the  incident  wave,  as  illustrated  in  Fig.  4.18 , could  be  further  exploited. 

If  a perfect  circularly  polarized  source  is  assumed  at  the  base  of  the 
ionosphere , the  polarization  of  the  received  signal  can  be  used  to  indicate 
the  elevation  angle  to  the  source;  the  more  nearly  circular  the  polarization, 
the  higher  the  elevation  angle . 

Another  mode  of  operation  has  been  tested , in  which  the  signals  from 
the  three  antenna  elements  are  recorded  on  a three-channel  tape  recor  der 
for  later  playback  arid  processing  by  the  frequency-tracking  direction-finder. 
This  method  requires  excellent  phase  and  amplitude  response  matching 
between  the  recorder  channels . However,  the  advantage  is  that  only  one 
direction-finding  instrument  is  required  to  analyze  data  from  any  number  of 
stations . Furthermore , more  detailed  analysis  of  signals  is  possible  by 
replaying  the  tape  witli  different  adjustments  of  the  instrument . Preliminary 
tests  have  shown  this  to  be  a feasible  approach . The  direction-finding 
errors  introduced  are  on  the  order  of  a few  degrees  for  a properly 


operating  recorder;  for  an  improperly  aligned  macidne , the  errors  may 
reach  5 to  15  degrees . 

Apart  from  the  direction-finding  function  the  frequency-tracking 
receiver  has  demonstrated  usehil  cape^ilities  in  the  unattended  detection 
and  analysis  of  VLF  activity . Unattended  reception  facilities  at  remote 
locations  or  on  satellites  employing  the  frequency-tracking  principle  can 
benefit  from  the  reduction  in  data  bandwidth  made  possible  by  automatic 
resd-time  extraction  of  frequency , amplitude , and  signal-to-noise  ratio 
information . 

A more  immediate  application  of  the  receiver  is  in  the  measurement  of 
the  growth  rate  of  the  artificially  stimulated  emissions  from  the  Siple  signal . 
Use  of  ordinary  filtering  techniques  is  hampered  by  the  changing  frequency 
of  the  emissions;  using  the  tracking  receiver,  a continuous  record  of  ttie 
amplitude  of  the  emission  can  be  produced  with  other  noise  filtered  out . 
Preliminary  experiments  along  this  line  have  been  carried  out  with  success. 
A minor  redesign  of  the  IF  filters  to  reduce  their  bandwidth  to  100  Hz  would 
optimize  the  instrument  for  this  application;  although  the  maximum  tracking 
rate  capability  would  be  reduced  by  a factor  of  ten , it  would  still  be  adequate 
for  the  relatively  slow  frequency  variations  of  the  Siple  signals . 
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APPENDIX  A . COMPARISON  OF  SMALL  AND  LARGE  APERTURE 

r-i~n — ^1 -1 — r i-ri--ri~-r iiifi  in -t 

DIRECTION-FINDERS 

Assume  two  model  direction-finding  sysf  - >ns , one  of  the  large-aperture 
interferometric  type  and  the  other  of  the  small-aperture  ^pe . 

The  interferometric  direction-finder  c^^nsists  of  two  pairs  of  antennas , 
each  pair  forming  a baseline  of  length  s . The  two  baselines  are  orthogonal . 
The  phase  difference  <p  between  the  two  antennas  constituting  each  baseline 
is  measured.  If  baseline  1 is  oriented  aiong  the  north-south  axis , and 
baseline  2 is  oriented  east-west , then 

= (2irA)8  cos  6 

<p^  = (2tt/,\,s  sin  & 

The  small  aperture  system  consists  of  crossed  loop  antennas , The 
amplitudes  of  the  signc  ls  from  the  two  elements  are  given  by 

. = A cos  8 
A„  = A sin  0 

t4 

In  the  presence  of  random  noise  and  measurement  errors,  the  effect  f 
on  t phase  measurement  (in  r&dians)  is  comparable  to  the  fractional  effect 
on  an  amplitude  measurement . That  is , 

(radians)  = AA/A  = e 


3i 


The  noisy  measurementa  produced  by  the  two  systeioa  are  therefore  ^ven 


(j>^  = (27rA)  s cos  6 + 


A.  = A ooe  8 + c,  A 
1 1 


4>^  = (2;rA)  8 sin  8 + 


X =>  A sin  8 + € A . 

2 A 


The  resolution  of  the  angle  of  arrival  from  the  sine  and  cosine  components 
measured  is  performed  as 


A /a  A 

8 = /a^  + jAg  . 


Substituting  the  above  quantities  and  rearranging, 


e = e' 


foom  which.it  can  be  seen  that  the  influence  of  the  noise  will  be  comparable 
when  (2ir/\)  a = 1,  or  s = \/2it . 


APPENDIX  B.  CIRCUIT  DIGRAMS  OF  THE  APPARATUS 

Schematic  diagrams  of  the  circuitry  used  to  implement  the  proto- 
type frequency-tracking  direction-finder  are  presented  herein . Fig . B . 1 
gives  a block  diagram  of  the  apparatus  showing  the  functions  performed 
on  each  circuit  card  and  the  interrelationships  between  the  cards.  The 
actual  schematics  of  the  cards , PCI  through  25  and  PC27 , are  given  in 
Figs . B . 2 through  B . 26 . The  backplane  wiring  of  the  card  cages  holding 
the  circuit  modules  is  diagrammed  in  B .27  and  B .28.  The  wiring  of  the 
control  panels  for  the  equipment  is  described  by  Figs . B . 29  and  B . 30 . A 
further  insight  into  the  operation  of  the  equipment  is  provided  by  the  drawings 
of  the  actual  front  panel  controls,  with  explanations,  in  Figs.  B.  31  and  B.32. 


FIGURE  B . I BLOCK  DIAGRAM  OF  THE  INSTRUMENT  SHOWING  CIRCUIT  CARD  FUNCTIONS 


FIGURE  B . 2 SCHEMATIC  OF  MASTER  LO  CARD , PCI . 


SCHEMATIC  OF  MIXER-VERT  CARD,  PC2. 


FIGURE  B .4  SCHEMATIC  OF  IP  FILTER- VERT  CARD , PCS 


FIGURE  B .5  SCHEMATIC  OF  UMITER  AMP  CARD , PC4. 


FIGURE  B . 6 SCHEMATIC  OF  DISCRIMINATOR  CARD  , PC  5 . 


FIGURE  B . 8 SCHEMATIC  OF  IF  FILTER  NOISE  CARD , PCT 


FIGURE  B . B SCtlEMATIC 


FIGURE  B . 10  SCHEMATIC  OF  LOG  nETEnvon  A irar^A 


FIGURE  B . 11  SCHEMATIC  OF  TRACK/RESE1’  LOGIC  CARD , PC  10 


FIGURE  B . 12  SCHEMATIC  OF  POTENTIOMETER  BUFFER  CARD , PCll 


PlQtJBE  B . 13  SCHEMATIC  OF  FREQ  TRAI^SLA.TOR  CARD . PC  12 


FIGORE  B.  14  SCHEMATIC  OF  NOISE  LO/MIXER  CARD,  PC  13. 


FIGURE  B . 15  SCHEMATIC  OF  PREAMP  CARD , PC  14 


FIGURE  B . 16  SCHEMATIC  OP  MIXER-NS  AND  EVI  CARD , PC  15 


FIGURE  B . 18  SCHEMATIC  OF  AGC  AMP  CARD . PC  18 


FIGUBJE  B . 19  SCHEMATIC  OF  SYNCHRONOUS  DETECTOR  CARD  , PC19. 


FIGURE  B . 20  SCHEMATIC  OF  DETECTOR  DRIVER  CARD  , PC20 


EW IWV  INPUTS 


FIGURE  B . 2 1 S CHEMATIC 


FIGU31E  B > 22  SCHEMATIC  OF  PHASE  METER  CARD  , PC22 


nOURE  B.23  SCHEMATIC  OF  FiVl  OSCILLATORS  CARD»  PC23. 


FIGURE  B . 24  SCHEMATIC  OF  EVENT  LOGIC  CARD  , PC24, 


FIGURE  B . 25  SCHEMATIC  OF  AUDIO  SUMMER  CARD » PC25 


FIGURE  B . 26  SCHEMATIC  OF  FREQUENCY  COUNTER  CARD , PC27 . 


m»-mou»  «oacETs  — «£a*  view 
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FIGUEE  B . 2S  SCHEMATIC  OP  TRACKING  RECEIVER  FRONT  PAlffiL , 


URE  B . 31  PICTOfUAL  DIAGRAM  OF  TRACKING  RECEIVER  FRONT  PANEL . 


Fimm  B , 32  PICTOBIAL  DIAGRAM  OP  DF  RECEIVER  FRONT  PAJJEL 
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NAA,  N?6,  GBR,  and  the  Omega  stations. & Well  defined  bearing  indications  were 
also  obtained  for  spherics  at  frequencf|s  in  the  1 tc  10  kilohertz  region. 
Successful  frequency-tracking  and  direc^on-finding  were  demonstrated  for  a 
wide  variety  of  whistler-mode  signals,  including  whistlers,  chorus,  and  dis- 
crete emissions.  The  technique  was  also  successfully  applied  to  emissions  in 
the  2 to  7 kilohertz  range  stimulated  by  the  VLF  transmitter  at  Siple  Station, 
Antarctica.  Direction  of  arrival  indications  on  whistlers  were  consistent 
with  the  duct  positions  predicted  by  frequency  profile  analysis.  Evidence 
was  also  presented  of  temporal  changes  in  the  direction  of  arrival  of  signals 
which  may  be  an  indication  of  duct  drifts.  This  instrument  provides,  for  the 
first  time,  a crntiruous,  imnediately  readable  record  of  the  direction  of 
arrival  of  whistler-mcde  signals.  It  could  also  be  applied  to  the  automatic 
acquisition  of  whistler-mode  signals  and  measurement  of  their  ajnplitude  and 
frequency  variations. 
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